AD-A021  924 


PROPAGATION  OF  A CHEMICAL  REACTION  THROUGH 
HETEROGENEOUS  LITH IUM-POLYTETRAFLUOROETHYLENE 
MIXTURES 

Richard  B.  Smith 
Pennsylvania  State  University 


Prepared  for: 

Naval  Sea  Systems  Command 


11  December  1975 


J 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 


Best  Avail 


3D’0 


mom 


0821  20 


; 

GNf 

05 


PROPAGATION  OF  A CHEMICAL  REACTION  THROUGH 
HETEROGENEOUS  LITHIUM-POLYTETRAFLUOROETHYLENE 
MIXTURES 


Richard  B.  Smith 


Oi 

o 


Technical  Memorandum 
File  No.  TM  75-302 
December  11,  1975 
Contract  No.  N00017-73-C-1418 


Copy  No.  (,j 


The  Pennsylvania  State  University 
Institute  for  Science  and  Engineering 
APPLIED  RESEARCH  LABORATORY 
Post  Office  Box  30 
State  College,  PA  16801 


REPRODUCED  BY 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 


U.S.  DEPARTMENT  OF  COMMERCE 
SPRINGFIELD,  VA.  22161 


c 


ju! 


J ^ly,1 


NAVY  DEPARTMENT 

NAVAL  SEA  SYSTEMS  COMMAND 


approved  For:  public  p:lease 
distribution  unlimited 


Best  Available  Co 


P) 


security  CLASSIFICATION  OF  THIS  PACE  r**Ti»n  Data  Enf»r»dl 


REPORT  DOCUMENTATION  PAGE 


. REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

TM  75-302 


4.  TITLE  fend  Sublllf) 

PROPAGATION  OF  A CHEMICAL  REACTION  THROUGH 
HETEROGENEOUS  LITH IUM-P  OLYTETRAFLUOROETHYLENE 
MIXTURES 


7.  AUTHORfJj 

Richard  B.  Smith 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT’S  CAT  ALOG  NUMBER 


S.  TYPE  OF  REPORT  A PERIOD  COVERED 

Ph.D.  Thesis,  March  1976 


6.  PERFORMING  ORG.  REPORT  NUMBER 

TM  75-302 


8.  CONTRACT  OR  GRANT  NUMBERftJ 

N00017-73-C-1418 


9.  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

The  Pennsylvania  State  University 

Applied  Research  Laboratory 

P.  0.  Box  30,  State  College,  PA  16801 


II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Naval  Sea  Systems  Command 
Department  of  the  Navy 

62 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  ft  WORK  UNIT  NUMBERS 


12.  REPORT  DATE 

December  11,  1975 


• 3.  NUMBER  OF  PAGES 


4.  MONITORING  AGENCY  NAME  ft  AOORESSf II  different  from  Confrotllna  Office) 

IS.  SECURITY  CLASS,  (of  thle  repo rt) 

Unclassified,  Unlimited 

15a.  OECL  ASSl  FI  CATION/ DOWNGRADING 
SCHEDULE 

16.  .DISTRIBUTION  STATEMENT  fof  (hit  Reporl) 


Approved  for  public  release,  distribution  unlimited, 
per  NSSC  (Naval  Sea  Systems  Command),  1/29/76 


17.  DISTRIBUTION  ST  ATEMENT  (of  the  ebutract  en  fared  In  Block  50,  If  different  from  Report) 


19-  KEY  WORDS  (Conf.'nu*  on  reverse  aide  if  neceasary  and  Identify  by  block  number) 


CHEMICAL  REACTION 
HETEROGENEOUS  MIXTURES 
LITHIUM  POLY TETRAFLU0R0 ETHYLENE 
PROPAGATION  VELOCITY 


REACTION  PREDICTION 
REACTION  RATE 
THERMAL  DECOMPOSITION 


1 20.  ABSTRACT  (Continue  on  reverse  aide  ii  necessary  andJdentify  by  block  number)  * 

An  analytical  and  experimental  investigation  of  the  propagation  of  a 
chemical  reaction  through  heterogeneous  lithium-poly tetraf luoroethylene 
(PTFE)  mixtures  was  conducted.  PTFE  in  the  form  of  both  solid  rod  and 
resin  was  reacted  with  lithium,  these  reactants  being  initially  at  room 
temperature  and  under  a vacuum.  The  reaction  sample  was  of  cylindrical 
geometry  with  the  PTFE  located  along  the  longitudinal  centerline.  The 
chemical  reaction  was  initiated  at  one  end  of  the  sample,  and  the  motion 


DD  1 JAN  73  1473  EDITION  OF  I NOV  65  15  OBSOLETE 


Best  Availabte1© 


UNCLASSIFIED 

ASSl  r IC  A T ION  OF  THIS  PAGE  (’H'hrn  fCnrrrcr* 


I 


UNCLASSIFIED 


j 

i 


J 

I 

4 


i 


4 

j 


j 


j 


I 


\ 


i 

( 


UNCLASSIFIED 


jccuhity  cl  rssinc  »tio«  or  this  * nctm™  n»i»  rnrnij 


ii 


ACKNOWLEDGMENTS 

The  author  wishes  to  express  his  sincere  appreciation  to 
his  advisor,  Professor  Donald  R.  Olson,  for  his  valuable  advice, 
timely  gu 'dance,  and  continued  encouragement  throughout  the  course 
of  this  investigation. 

he  also  desires  to  thank  Dr.  Laurence  A.  Ferris,  of 
E.  I.  duPont  de  Nemours  and  Company,  for  his  kind  assistance  in 
obtaining  a portion  of  the  PTFE  used  in  the  experimental  program. 

The  author  wishes  to  acknowledge  the  support  for  this 
investigation  by  the  Applied  Research  Laboratory  of  The  Pennsylvania 
State  University  under  contract  with  the  Naval  Ordnance  Systems  Command. 


V 


iii 


TABLE  OF  CONTENTS 


ACKNOWLEDGMENTS 


LIST  OF  TABLES v 

LIST  OF  FIGURES vi 

NOMENCLATURE.  vii 

I.  INTRODUCTION  1 

1.1  General  Statement  of  the  Problem 1 

1.2  Some  Similar  Reaction  Studies  Involving  One  or 

More  Condensed  Phases  9 

1.2.1  Lithium-Gas  Surface  Reactions 10 

1.2.2  Composite  Solid  Propellant  Combustion.  ...  13 

1.2.3  Flame  Spread  Over  the  Surface  of  Condensed 

Phase  Materials.  ......  14 

1.2.4  Exothermic  Alloying  of  Metals 17 

1.2.5  Combustion  of  Magnesium-PTFE  Pyrotechnic 

Pellets 18 

1.3  Specific  Research  Objectives 19 

II.  THEORETICAL  CONSIDERATIONS  21 

2.1  Description  of  the  Proposed  Reaction  Mechanism.  . . 21 

2.2  One-Dimensional  Model  22 

2.3  Two-Dimensional  Model  23 

2-4  Details  of  the  Numerical  Solution 31 

2.5  Results  of  the  Analytical  Model  40 

III.  EXPERIMENTAL  INVESTIGATIONS 51 

3.1  Experimental  Apparatus 51 

3.2  Experimental  Procedure 59 

IV.  EXPERIMENTAL  AND  THEORETICAL  RESULTS  65 

4.1  Experimental  Conditions  ...  65 

4.2  Test  Runs  with  9.52  mm  Diameter  PTFE  Samples.  ...  67 

4.3  Test  Runs  with  6.35  mm  Diameter  PTFE  Samples.  ...  76 

4.4  Test  Runs  with  PTFE  Samples  of  Other  Diameters.  . . 85 

4.5  Experimental-Theoretical  Comparisons 89 

V.  SUMMARY 96 


TABLE  OF  CONTENTS  (CONTINUED) 


BIBLIOGRAPHY 99 

APPENDIX  A:  DERIVATION  OF  THE  FINITE  DIFFERENCE 

EXPRESSIONS 103 

APPENDIX  B:  COMPUTER  PROGRAM  FLOW  DIAGRAM 113 

APPENDIX  C:  VALUES  OF  PHYSICAL  AND  CHEMICAL  PROPERTIES 

USED  IN  THE  MODEL 116 


V 


LIST  OF  TABLES 


Table 

Title 

Page 

1 

Summary  of  Lithium-PTFE  Test 
Diameter  PTFE  Resin  Samples. 

Runs 
• * 

with 
« • • 

9.52  mm 

. . 68 

2 

Summary  of  Lithium-PTFE  Teat 
Daimeter  PTFE  Resin  Samples. 

Runs 

with 

t » » 

6.35  mm 

. . 77 

3 

Summary  of  Lithium-PTFE  Test 
Diameter  PTFE  Rod  Samples.  . 

Runs 

with 
• * • 

6 . 35  mm 

. . 78 

<♦ 

Summary  of  Lithium-PTFE  Test 
PTFE  Rod  Samples  ...... 

Runs 

with 

» • i 

Commercial 

. . 86 

5 

Summary  of  Results  of  Successful  ' 
Test  Runs  

Lithium-PTFE 

. 92 

6 

Values  of  Physical  and  Chemical  Properties  Used 
in  the  Model 

. 116 

vi 


LIST  OF  FIGURES 

Figure  Caption  Page 

1 Preliminary  Test  Apparatus.  ...  5 

2 Improved  Test  Apparatus  « ........  7 

3 Lithium-PTFE  Reaction  Model 24 

4 Sketch  of  Finite  Difference  Grid  Network 33 

5 V*  as  a Function  of  X*  for  Two  Values  of  L*/H* 

for  PTFE-Lithium  Mass  Ratio  » 0.379  43 

6 Final  Equilibrium  Temperature  as  a Function  of 

PTFE-Lithium  Maas  Ratio  .......  45 

7 IX  Tolerance  Dimensionless  Steady  Reaction 
Propagation  Velocity  as  a Function  of  PTFE-Lithium 

Mass  Ratio.  47 

8 Dimensionless  Temperature  Profiles  at  2.5X  Tolerance 
Steady  Conditions  for  L*/H*  * 1.7  and  PTFE-Lithium 

Mass  Ratio  • 0,206.  . 48 

9 Sketch  of  Test  Apparatus 52 

10  Photograph  of  Test  Apparatus, 54 

11  Sketch  of  Test  Facility  ...  . , . . 55 

12  Photograph  cf  Test  Facility  ......  57 

13  Photcgraph  of  Sample  Tube,  Thermocouple  Probes,  and 

Core  Rod 60 

14  Sketch  of  Temperature-Time  Record  for  a Typical 

Test  Run-  . 62 

15  Sketch  of  Finite  Difference  Grid  Patterns 


. 104 


~+y  ■ •^r*nsr*T nmi- 


vii 


NOMENCLATURE 


Description 


Dimensionless  parameter  defined  by  Equation  (2.15) 

Frequency  factor  in  PTFE  thermal  decomposition  rate  law. 
Equation  (2.4)  fs“^) 

Dimensionless  constants  defined  by  Equations  (A. 7),  (A. 17), 
(A. 25)  and  (A. 29) 

Specific  heat  at  constant  pressure  (J/kg  - K) 

Dimensionless  parameter  defined  by  Equation  (2.15) 

Activation  energy  in  PTFE  thermal  decomposition  rate  law, 
Equation  (2.4)  (J'mol) 

Natural  logarithm  base 

Dimensionless  parameter  defined  by  Equation  (2.15) 
Dimensionless  thickness  cf  lithium-PTFE  reaction  sample 
Dimensionless  thickness  of  PTFE 
Standard  heat  of  reaction  (J/kg  reactants) 

Thickness  (m) 

Constant  defined  by  Equation  (4.1)  (s) 

Thermal  conductivity  (W/m-K) 

Dimensionless  length  ct  lithium-PTFE  reaction  sample 
Length  of  lithium-PTFE  reaction  sample  (m) 

Dimensionless  PTFE-lithium  mass  ratio 
Mass  of  PTFE  (kg) 

Initial  mass  of  PTFE  (kg) 

Number  of  Y*  direction  grid  points  in  PTFE 

Number  of  Y*  direction  grid  points  in  stoichiometric  lithium 


Number  of  X*  direction  grid  points  in  lithium-PTFE  reaction 
sample 


NOMENCLATURE  (CONTINUED) 


Description 

3 

Rate  of  heat  generation  per  unit  volume  (W/m  ) 

Dimensionless  parameter  defined  by  Equacion  (2.15) 

Universal  gas  constant  (J/mol  - K) 

Dimensionless  temperature  defined  by  Equation  (2.11) 
Dimensionless  initial  temperature 
Temperature  (K) 

Adiabatic  flame  temperature  of  lithium-PTFE  reaction  at  one 
atmosphere  total  pressure  (lO 

Initial  temperature  (K) 

Time  (s) 

Characteristic  time  required  to  tnermally  decompose  90X 

of  an  initial  amount  of  PTFE  at  T (s) 

at 

Dimensionless  reaction  propagation  velocity  defined  by 
Equation  (2.26) 

Characteristic  velocity  defined  by  Equation  (2.27)  (m/s) 

Ratio  cf  mass  of  reactants  consumed  tc  mass  of  PTFE  thermally 
decomposed 

Dimensionless  coordinate  parallel  with  direction  of  reaction 
front  propagation,  defined  by  Equation  (211) 

Dimensionless  X*  direction  grid  spacing 

Number  of  moles  of  lithium  present  in  reaction  sample  in 
excess  of  that  required  by  stoichiometry 

Coordinate  parallel  with  direction  of  reaction  front 
propagation  (m) 

Dimensionless  coordinate  perpendicular  to  direction  of 
reaction  frcnc  propagation,  defined  by  Equation  (2.11) 

Dimensionless  Y*  direction  grid  spacing 

Coordinate  perpendicular  tc  direction  of  reaction  front 
propagation  (m) 


NOMENCLATURE  (CONTINUED) 


Symbol  Description 

Z*  Mass  fraction 

2 

a Thermal  diffueivity  (m  /s) 

3 

p Density  (kg/m  ) 

T*  Dimensionless  time  defined  by  Equation  (2.11) 


Subscripts 


i 

Pertains 

to 

a specific  grid  point 

J 

Pertains 

to 

a specific  grid  point 

L 

Pertains 

to 

lithium 

LS 

Pertains 

to 

stoichiometric  lithium 

LT 

Pertains 

to 

lithium-PTFE  interface 

LX 

Pertains 

to 

excess  lithium 

P 

Pertains 

to 

reaction  products 

S 

Pertains 

to 

steady  conditions 

T 

Pertains 

to 

PTFE 

Superscript 


N 


Pertains  to  dimensionless  time  level 


CHAPTER  I 


INTRODUCTION 

1 . 1 General  Statement  of  the  Problem 

The  past  decade  has  witnessed  a vast  intensification  in  the 
search  for  new  sources  of  energy  in  the  form  of  heat.  In  many 
applications  such  as  electric  power  generation,  supply,  cost,  and 
environmental  considerations  are  of  primary  concern  in  the  selection 
of  a heat  source.  Space,  underwater,  and  several  other  applications, 
however,  dictate  that  high  energy  density  is  of  the  utmost  importance, 
energy  density  being  defined  as  the  amount  of  heat  generated  per  unit 
weight  of  fuel  and  oxidizer  consumed.  The  chemical  reactions  between 
the  alkali  metals  and  the  halogen  gases  are  all  quite  attractive  from 
an  energy  density  standpoint,  and  of  these  several  possible  fuel- 
oxidizer  combinations,  the  reaction  between  lithium  and  fluorine 
produces  the  largest  amount  of  heat  per  unit  weight  of  reactants. 

Lithium  is  the  lightest  known  metal,  and  is  relatively 
plentiful.  Although  it  is  extremely  reactive  at  elevated  temperatures, 
lithium  may  be  handled  with  reasonable  ease  and  safety  at  temperatures 
up  to  slightly  above  its  melting  point,  454  K.  Fluorine,  on  the 
other  hand,  is  exceedingly  toxic  and  very  difficult  to  handle  under 
all  conditions.  The  highly  fluorinated  polymers  such  as 
polytetraf luoroethylene , PTFE  (Teflon*,  Fluon**,  Halon***),  however. 


*Pesistered  trademark,  E.I.  duPont  de  Nemours  and  Company 
**Registered  trademark,  ICI  America  Incorporated 
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provide  a convenient,  safe,  and  inexpensive  source  of  fluorine. 

PTFE,  chemical  formula  nC^F^,  is  76t  fluorine  by  weight,  and  in  its 
solid  form  is  slightly  more  dense  in  fluorine  than  is  liquid  fluorine 
at  its  normal  boiling  point.  Available  commercially  in  the  form  of 
rod,  sheet,  and  resin,  PTFE  is  known  primarily  for  its  chemical 
inertness.  Product  literature  indicates  that  PTFE  will  react  with  the 
molten  alkali  metals,  although  information  concerning  the  conditions 
under  which  the  reaction  might  be  initiated  and  the  order  of 
magnitude  of  the  attendant  rate  is  not  available.  Were  it  possible  to 
initiate  the  lithium-PTFE  reaction  without  undue  effort,  and  were  this 
reaction  to  proceed  at  a reasonably  rapid  rate,  this  reactant 
combination  might  well  prove  attractive  in  applications  requiring  a 
high  energy  density  heat  source.  Therefore,  some  preliminary  testing 
was  undertaken  in  order  that  the  answers  to  these  questions  might  be 
ascertained. 

At  the  outset  of  preliminary  testing,  it  was  decided  that 
the  configuration  in  which  the  lithium-PTFE  reaction  might  be  examined 
with  the  most  generality  was  that  of  a flat  circular  wafer  of  lithium 
covered  by  the  PTFE  oxidizer.  It  was  hoped  that  with  the  sample  under 
an  initial  vacuum,  the  reaction  could  be  ignited  at  the  center  of 
the  wafer  by  means  of  an  electrically  heated  wire,  and  allowed  to 
propagate  radially.  Accordingly,  some  feasibility  tests  were  run  using 
a 25.4  mm  diameter  lithium  wafer  covered  with  PTFE  resin.  Although 
a considerable  number  of  attempts  at  igniting  the  reaction  with  the 
electrically  heated  wire  alone  proved  fruitless,  it  was  found  that  by 
augmenting  the  hot  wire  with  a small  amount  of  monopropellant,  the 
reaction  could  be  ignited  quite  reliably.  The  volume  cf  the 
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containment  vessel  used  in  these  tests  was  only  slightly  greater  than 
that  required  by  the  reactants;  very  little  free  volume  was  present 
above  the  sample. 

The  use  of  this  experimental  configuration  having  appeared 
feasible,  a larger  scale  test  apparatus  to  accommodate  a 140  mm 
diameter  lithium  sample  was  constructed.  In  this  device,  a rather 
large  free  volume  above  the  reaction  sample  was  provided  in  order  that 
test  preparation  might  be  facilitated.  After  the  lithium  wafer  and  the 
PTFE  resin  were  installed,  a thin  stainless  steel  disc  was  placed  on 
top  of  the  PTFE  resin.  Thermocouples,  which  were  intended  to  be  used 
to  monitor  the  position  of  the  propagating  reaction  front,  were 
inserted  through  holes  placed  in  this  disc  at  various  radial  and 
circumferential  locations.  The  thermocouple  measuring  junctions  were 
located  at  the  lithium-PTFE  interface.  A larger  hole  was  located  at 
the  center  of  the  disc  to  provide  an  area  for  the  ignition  of  the 
reaction.  With  this  device,  all  attempts  at  initiating  a self- 
propagating  reaction  under  a vacuum  with  the  monopropellant-augmented 
hot  wire  proved  unproductive.  In  some  cases  the  PTFE  resin  was  seen 
to  be  blown  away  from  the  ignition  area  by  the  burning  of  the 
monopropellant;  in  others,  the  reaction  appeared  tc  have  propagated 
a short  distance  and  then  stopped.  In  all  cases,  a considerable  amount 
of  material  thought  to  1 e reaction  products  from  the  monopropellant 
was  found  in  the  free  volume  over  the  sample  after  ignition  was 
attempted . 

From  the  results  of  these  tests,  it  was  concluded  that  the 
lithium-PTFE  reaction  would  not  self-propagate  in  a radial  manner,  in 
which  the  surface  area  covered  by  the  reaction  front  was  forced  to 
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increase  during  the  propagation.  For  this  reason,  it  was  decided 
to  continue  preliminary  testing  in  a cylindrical  configuration,  in 
which  a constant  reaction  front  surface  ^'ca  was  maintained. 
Additionally,  it  was  felt  that  the  use  of  a commercially  available 
igniter  in  place  of  the  monopropellant-augmented  hot  wire  would 
substantially  Improve  the  probability  of  a reliable  ignition  of  the 
reaction , 

An  initial  experiment  involving  lithium  and  PTFE  resin  in  the 
cylindrical  configuration  was  conducted  in  an  apparatus  as  shown  in 
Figure  1.  The  device  consisted  of  a length  of  stainless  steel  tubing 
with  a welded  closure  at  one  end  and  a tube  fitting  at  the  other.  A 
commercially  available  lithium  rod  of  25.4  mm  diameter  and  203  mm 
length  was  inserted  into  this  tube,  and  a 9.5  m diameter  hole  along 
the  longitudinal  axis  of  the  lithium  rod  was  produced  by  driving  a 
length  of  tubing  through  the  lithium  and  removing  the  resulting 
slug.  A quantity  of  PTFE  resin  was  tamped  into  this  hole,  and  an 
igniter  was  positioned  at  the  upper  end  frr  the  purpose  of  initiating 
the  reaction  While  these  preparations  were  being  carried  out,  an 
attempt  was  made  to  continuously  bathe  the  sample  in  argon  gas  to 
minimize  air  contamination  of  the  lithium.  Thermocouples  welded  to 
the  outside  of  the  tube  along  Its  length  were  connected  to  a recording 
oscillograph  in  order  that  the  motion  of  the  reaction  front  might  oe 
monitored  Thermochemical  calculations,  based  on  the  assumption  of 
lithium  fluoride  and  lithium  carbide  reaction  products,  indicated 
that  sufficient  excess  lithium  was  present  to  limit  the  final  mean 
temperature  of  the  material  within  the  tube,  after  the  reaction  had 
taken  place,  to  a reasonable  level  of  1100-1150  K.  Before  the  test 


Figure  1 Preliminary  Test  Apparatus 
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was  run,  the  device  was  evacuated  and  the  valve  was  closed,  isolating 
the  system.  The  igniter  was  then  fired,  and  the  reaction  was  seen  to 
traverse  the  length  of  the  sample  in  approximately  0.035  seconds, 
corresponding  to  an  average  velocity  of  about  5.5  m/s.  The  arrival 
of  the  reaction  front  at  each  subsequent  thermocouple  location  was 
evidenced  by  a sharp  deflection  on  the  corresponding  analog  output 
record,  indicating  a rapid  increase  in  temperature  at  that  point. 
Visual  inspection  of  the  apparatus  and  its  ccntents  after  the  test 
run  indicated  that  all  of  the  PTFE  resin  had  been  consumed  in  the 
reaction.  In  addition,  no  evidence  of  high  pressure  within  the  device 
was  seen.  It  was,  however,  recognized  that  the  placement  of  the 
thermocouples  on  the  outside  of  the  sample  tube  introduced  a 
considerable  time  lag  in  the  measuring  system.  It  was,  therefore, 
decided  to  carry  out  some  more  refined  tests  in  an  attempt  to  get  some 
more  basic  information  on  the  characteristics  or  this  reaction. 

The  second  series  of  test  runs  in  the  cylindrical 
configuration  was  carried  out  in  an  apparatus  as  shown  in  Figure  2 
The  sample  tube  was  again  of  sufficient  size  to  allow  the  insertion 
of  a commercial  25.4  mm  diameter  x 203  mm  lithium  red,  and  the  hole 
into  which  the  PTFE  vesin  was  placed  was  again  created  by  forcing  a 
length  of  tubing  through  the  lithium  sample.  In  this  configuration, 
however,  the  igniter  was  threaded  into  a plug  which  was  in  turn 
weldeo  into  one  end  of  the  sampJe  tube,  and  the  closure  at  the 
opposite  end  of  the  tube  was  affected  by  a removable  cap  incorporating 
a redundant  high  temperature  0-ring  seal-  This  cap  was  held  on  during 
testing  by  a retainer  riug  and  four  tie  bolts.  The  most  important 
refinement  included  in  this  device  was,  however,  che  placement  of  the 


8 


thermocouple  measuring  junctions  at  the  lithium-PTFE  interlace.  This 
was  accomplished  through  the  use  of  fine  wire  '0  i'>  mm  diameter) 
exposed  junction  thermocouple  probes  having  a sheath  diameter  of  1.0 
mm.  The  thermocouple  probes  were  installed  after  the  lithium  rod  was 
Inserted  and  the  hole  for  the  PTFE  resin  was  created.  They  were 
sealed  at  the  surface  of  the  sample  tube  with  high  temperature 
epoxy.  After  the  thermocouple  probes  and  the  PTFE  resin  were  in 
place  and  the  cap  and  tie  bolts  secured,  the  thermocouple  probes 
were  connected  to  a recording  oscillograph,  and  the  device  was 
evacuated  in  preparation  for  the  firing  of  the  igniter.  During  the 
period  of  time  that  the  lithium  was  exposed  to  the  air,  the  device 
was  continuously  bathed  in  argon  gas  in  an  attempt  to  minimize 
contamination  of  the  specimen.  A total  cl  four  successful 
preliminary  tests  were  conducted  with  this  apparatus  The  nominal 
amounts  of  reactants  involved  were  50  g of  lithium  ar.d  12  g of  PTFE 
resin,  and  the  average  rate  of  propagation  or  the  reaction  rront  was 
36. 3 m/s  with  a range  of  29  to  <»5  m's.  On  additional  test  was 
conducted,  in  which  the  reaction  was  seen  to  propagate  o«e: 
approximately  10%  of  the  length  of  the  specimen  and  then  stop.  This 
was  attributed  to  non-uniform  packing  of  the  PTFE  re9in 

From  the  results  of  these  experiments  in  the  cylindrical 
configuration,  it  was  concluded  that  the  chemical  reaction  between 
lithium  and  PTFE  resin,  each  initially  at  rccm  temperature  and 
under  a vacuum,  was  self-piopagatmg  when  p-operLy  ignited,  and  was 
extraordinarily  fast  It  was  therefore  decided  that  a more  thorough 
investigation  into  this  reaction  process  should  be  undertaken. 
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although  It  was  clear  that  the  ignition  of  this  reaction  was  a 
difficult  task  requiring  an  ignition  source  high  in  both  quantity  and 
rate  of  energy  output 

1. 2 Some  Similar  Reaction  Studies  Involving  One  or  Moie 

Condensed  Phases 

The  chemical  reaction  between  lithium  and  PTFE,  which  is 
heterogeneous  on  both  the  macoscopic  and  the  microscopic  scale  in 
the  configuration  utilized  in  the  preliminary  expeciments  described, 
is  believed  to  follow  the  stoichiometry 

6Li  + C,F.  - 4LiF  + Li„C_  ' . 1 > 

2 4 2 2 

At  room  temperature,  the  reactants  are  in  the  solid  prase,  as  are  the 
products.  An  estimate  of  the  adiabatic  tlame  temperature  of 
stoichiometric  amounts  of  lithium  and  PTFE,  based  on  cuoienciy 
available  thermochemical  information  indicated  u.at  this 

temperature  was  of  the  order  of  2000  K,  well  abote  the  meeting  point 
cf  standard  materials  of  construction  su:h  as  stainless  steel  In 
order  to  preclude  the  melting  of  the  sample  ccncaine*  tube  during 
testing,  a quantity  of  lithium  in  ex  .ess  of  that  requi-ed  tc  combine 
stoi chiometrically  with  the  PTFE  present  was  in~l-.dcd  The  oveiai- 
reaction  may  thus  be  represented  by 

'6  * XS'  li  C.F  * •‘.LI*  + li~C.  + XS  11  , a.?) 

Z u z z 

where : 

XS  * number  of  mcles  of  lithium  poesent  in  excess  of  chat 


required  by  stoichiometry 


10 


In  the  preliminary  experiments,  the  amount  of  excess  lithium  was 
chosen  such  that  the  mean  temperature  of  the  reaction  products  and 
excess  lithium  immediately  after  the  reaction  had  taken  place  was 
1100*1150  K.  The  melting  temperature  of  lithium  fluoride  is  1121  K 
[1];  lithium  carbide  does  not  melt  but  decomposes  to  graphite  and 
lithium  vapor  at  1 atmosphere  pressure  at  182S  K [2).  The  material 
within  the  sample  container  tube  immediately  after  the  reaction  had 
occurred  thus  consisted  of  liquid  lithium,  solid  or  liquid  LiF,  and 
solid  Li2C2. 

A literature  search  was  conducted  on  the  subject  of 
macroscopically  heterogeneous  chemical  reactions  in  which  condensed 
phase  reactants  form  condensed  phase  products,  from  which  no 
pertinent  references  were  obtained.  It  was  therefore  decided  that  an 
examination  of  some  previous  reaction  studies  involving  one  or  more 
condensed  phases  might  prove  instructive. 

1.2.1  Lithium-Gas  Surface  Reactions 

The  chemical  reaction  between  dry  cxygen  and  a molten  lithium 
surface  was  Investigated  by  Tyzack  and  Lcr.gton  [4]  at  subatmospher l c 
pressures  and  at  temperatures  up  to  1173  K The  primary  purpose  of 
this  work  was  to  measure  the  ignition  cemperature  of  the  lithium- 
oxygen  reaction  as  a function  cf  oxygen  pressure.  The  minimum  ignition 
temperature  observed  was  900-920  K which  oczurred  at  pressures  from 
0.13  to  0,67  bar.  Ignition  temperatures  as  high  as  1073  K were  found 
at  pressures  above  and  below  this  range  Since  the  reaction  product, 
lithium  oxide,  was  in  the  solid  phase  at  the  temperatures  encountered 
in  t e experiments,  some  reaction  rate  data  was  obtained  by  measuring 
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Che  change  in  the  oxygen  reservoir  pressure  as  a i unction  of  time  over 
a period  of  approximately  30  minutes  following  ignition  Only 
representative  rate  data  were  reported,  and  no  estimates  of  activation 
energy  or  rate  constant  were  presented.  No  analytical  model  was 
proposed 

Longton  (5]  studied  the  surface  reaction  between  molten 
lithium  and  nitrogen  over  the  pressure  range  0 07-1.0  bar,  and  at 
temperatures  from  57^  to  723  K.  From  experiments  similar  to  those  of 
Tytack  and  Longton  (4],  but  with  a typical  time  duration  or  100 
minutes,  this  reaction  was  iound  to  follow  a parabolic  rate  law  after 
a short  initial  period.  The  activation  energy  was  determined,  *nd 
the  rate  constant  was  tabulated  as  a function  ct  lithium  temperature 
and  nitrogen  pressure.  No  analytical  effort  was  reported. 

Chandrasekharaiah  and  Margrave  [6]  conducted  investigations 
of  the  oxidation  and  nltridation  of  a lithium  surface  In  this  study, 
no  oxidation  race  measurements  were  made  sln<e  the  rates  were  found  to 
be  immeasurably  small  at  row  temperatures,  and  no  containment  material 
could  be  found  which  was  suitable  for  use  at  high  temperature 
Nltridation  rate  measurements  were  made  during  30  to  j.  20  minute  time 
periods  following  exposure  of  the  metal  tc  nitrogen  gas  at  pressures 
of  0.15  and  027  bar  The  metal  surface  temperatures  du'ing  these 
experiments  ranged  from  512  to  583  K.  The  measu'ed  reaction  rates  were 
erratic  and  the  apparent  activation  energy  was  reported  to  lie  between 
42  and  188  kj'mole.  Nc  analytical  mode*  was  postulated 

i lie  nltridation  ci  solid  and  mclten  lithium  was  studied  by 
McFarlane  and  Tompkins  (71  at  temperatures  from  273-638  K and  at 
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pressures  from  0.13-0.67  bar.  The  reaction  rate  of  nitrogen  with  a 
solid  lithium  surface  was  observed  to  increase  with  time  until  the 
point  at  which  a critical  thickness  ot  the  lithium  nitride  reaction 
product  had  formed  Subsequent  to  the  attainment  of  this  critical 
thickness,  the  reaction  rate  was  reported  to  be  constant  with  time  but 
pressure  dependent  A pressure-independent  parabolic  rate  with 
64.9  kJ/tnole  activation  energy  was  observed  tor  the  surface  reaction 
of  nitrogen  with  molten  lithium.  Test  data  were  obtained  during 
periods  of  up  to  200  minutes  following  exposure  of  the  lithium  te  the 
nitrogen,  and  no  analytical  model  was  presented 

The  surface  reaction  of  molten  lithium  with  nitrogen  and 
oxygen,  as  well  as  with  the  halogenated  gases  CC1_F  C . F0  and  SF,  was 
investigated  by  Little  [8].  Rate  data  were  presented  tor  the  reaction 
of  lithium  with  001^2  and  C^Fg  at  times  trom  0,5  tc  300  seconds 
following  exposure  of  the  clean  metal  surface,  at  temperatures  from 
640  to  1000  K,  to  the  reactant  gas  at  pressures  f'om  0.001  to  0.1  bar. 
The  reaction  rate  between  lithium  and  nitrogen  was  'epr-ced  to  be 
pressure-independent  while  the  rates  of  the  tea:  ions  ct  lithium  with 
both  CC1  F and  C F were  observed  to  depena  cn  the  gas  pressure 
The  ignition  tempecatu’e  of  the  lithium  oxygen  reaction  was  reported 
tc  be  880  K,  and  the  lithium-SF^  reaccicn  was  observed  to  ignite  at 
1063  K,  corresponding  tc  the  temperature  a*,  which  the  reaction  produces 
become  liquid.  A criterion  fer  the  occurrence  of  a surface  reaction 
was  developed,  but  no  analytical  reaction  rate  model  was  presented 
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1.2.2  Composite  S-lid  Propellant  Combustion 

Composite  solid  propellants  are  plastic-like  materials 
consisting  of  small  oxidizer  particles  embedded  in  a fuel  matrix. 
Ammonium  perchlorate  is  the  most  commonly  used  oxidizer,  and  the  fuel, 
which  also  acts  as  a binder,  is  generally  a plastic  or  rubber  material. 
Composite  solid  propellants  are  therefore  homogeneous  on  the 
macroscopic  scale,  but  microscopically  heterogeneous.  Typical 
experimental  investigations  (9}  involve  the  burning  of  the  propellant 
at  constant  pressure  over  the  range  7-200  bar ■ The  products  of 
combustion  are  gaseous.  Linear  self-propagating  burning  rates  of  the 
order  of  10-100  ram/s  are  found  to  be  a function  of  propellant 
composition  and  particle  size,  environmental  pressure,  propellant 
temperature,  and  the  nature  of  the  gas  flow  relative  tc  the  burning 
surface.  The  actual  combustion  process  is  very  complex,  with 
endothermic  decomposition  reactions  supplying  gaseous  fuel  and 
oxidizer  which  subsequently  participate  In  one  or  more  exothermic 
chemical  reactions. 

The  earliest  analytical  treatment  of  composite  solid 
propellant  combustion  is  the  "Granular  Diffusion  Flame"  model  :t 
Summerfield,  et  al.,  (10).  Through  a primarily  qualitative  line  of 
reasoning,  the  authors  were  able  to  produce  a relationship  between 
the  linear  burning  rate  and  the  environmental  pressure.  Fenn's 
"Phalanx  Flame"  model  [11],  which  is  also  qualitative  in  nature, 
predicts  a similar  burning  rate-pressure  relationship  throu  gh  a 
different  set  cf  arguments.  Both  the  Summerfield,  et  al . , (10), 
and  the  Fenn  (11)  models  are  essentially  heuristic  devices,  intended 
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Co  provide  some  physical  Insight  into  the  phenomena  occurring,  and 
to  serve  as  a starting  point  for  the  development  of  more  complete 
theories,  .lore  rigorous  analytical  treatments  by  Hermance  (12], 
and  Culick  and  Dehority  [13]  have  also  been  reported,  but  as  yet 
no  widely  accepted  comprehensive  theory  of  the  combustion  of 
composite  solid  propellants  exists, 

1.2.3  Flame  Spread  Over  the  Surface  of  Condensed  Phase  Materials 

An  understanding  of  the  mechanism  by  which  a flame  spreads 
over  the  surface  of  a solid  or  liquid  material  in  an  oxidizing 
gaseous  environment  is  of  major  importance  in  such  areas  cf  interest 
as  the  ignition  of  solid  propellant  rocket  motors  and  the  advance  of 
forest  fires.  In  the  case  in  which  the  surface  of  a soild  is 
involved,  the  'lame  spreading  velocity  has  been  found  experimentally 
[14-20]  to  be  effected  by  the  chemical  composition  and  surface 
texture  of  the  condensed-phase  material,  the  details  ci  the  movement 
(if  any)  of  the  gaseous  environment,  the  orientation  and  geometry 
of  the  surface,  t*'e  thickness  cf  the  fuel  mass,  the  ambler:  oxidizer 
concentration,  and  the  total  ambient  pressure.  Analytical  models  of 
cne  spread  of  a flu.ie  over  a solid  surface  differ  chiefly  in  che 
assumptions  made  relative  to  che  dominant  mode  of  heat  transfer  iron 
the  advancing  flame  forward  to  preheat  and  vaporize  che  unbcrned 
surf  ace . 

McAlevy  and  Magee  [21]  postulated  that  the  flame  spreading 
velocity  is  controlled  by  the  gas-phase  heat  and  mass  transfer 
occurring  within  a small  "ignition  region"  immediately  ahead  of  the 
flare  Through  a simplitied  analytical  treatment  and  the  use  of 
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experimentally-determined  temperature  profiles,  the  authors  arrived 
at  a power-law  relationship  between  the  total  pressure,  the  oxidizer 
mole  fraction,  and  the  flame  spreading  velocity  which  was  seen  to 
correlate  experimental  data  with  a reasonable  degree  of  success. 
Solid-phase  phenomena  were  ignored.  In  a late''  effort  by  Lastrina, 
Magee,  and  McAlevy  [22],  both  gas  and  solid-phase  analyses  were 
included.  Flame  spreading  velocity  control  by  the  processes  taking 
place  within  the  ignition  region  was  again  postulated,  but 
differences  in  these  controlling  processes  caused  by  the  fuel  being 
thermally  "thick"  (temperature  gradients  within  the  fuel  in  the 
direction  normal  to  the  surface  are  non-negligible ) , or  thermally 
"thin"  (temperature  of  the  fuel  is  constant  in  the  direction  normal 
to  the  surface)  were  taken  into  account.  Simple  algebraic  equations 
of  the  power-law  type  for  the  flame  spreading  velocity  over  thermally 
thick  and  thermally  thin  fuel  beds  resulted.  These  relations, 
although  requiring  experimental  information  in  order  to  be  evaluated, 
have  been  used  quite  successfully  by  the  authors  to  correlate 
experimental  results  - 

Flame  spreading  over  both  thermally  thick  and  thermally 
thin  solid  fuel  beds  has  also  been  treated  analytically  by  de  Rls 
[23],  It  was  assumed  that  the  heat  transfer  from  the  flame  forward 
to  the  unburned  fuel  was  by  solid  and  gas-phase  conduction  in  the 
case  of  the  thermally  thick-fuel,  and  by  gas-phase  conduction  only 
in  the  thermally  thin  situation.  By  the  choice  not  tc  identify  a 
finite  ignition  region,  and  the  assumption  of  infinite  reaction 
rates,  the  model  was  made  amenable  to  treatment  by  ditiusion  flame 
theory,  and  an  exact  mathematical  solution  was  possible.  The 
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resulting  expressions  for  the  flame  spreading  velocity,  vuich  also 
require  the  input  of  experimental  data  in  order  to  be  evaluated,  ha^e 
been  shown  to  be  similar  [15]  to  those  developed  by  Lastrina, 

Magee,  and  McAlevy  (22)  for  fuel  beds  of  like  thermal  thickness. 

This  similarity  is  surprising  in  light  of  the  considerable  difference 
in  assumptions  regarding  the  existence  f the  finite  ignition  region 

Another  model  of  the  process  of  flame  spreading  over  solid 
surfaces,  in  which  radiation  is  assumed  to  be  the  dominant  mode  of 
heat  transfer  from  tne  flame  to  the  unburned  fuel,  has  been  proposed 
by  Tarifa,  Notario,  and  Torralbo  [24]  By  ignoring  conduction  and 
convection  of  energy,  a consideration  of  the  effects  of  the  thermal 
thickness  of  the  fuel  bed  was  made  unnecessary.  The  resulting 
implicit  expressions  involving  the  flame  spreading  velocity  are 
sufficiently  complex  uS  to  preclude  quantitative  comparison  with 
experimental  results. 

Heat  conduction  forward  through  the  solid  fuel  bed  is 
assumed  to  be  the  major  mode  of  heat  transfer  in  an  additional  mcdel 
proposed  by  Fernandez-Pello  and  Williams  [1')-  This  model  Vis 
specifically  constructed  for  the  spread  ci  a laminar  flame  over  the 
surface  of  polymethylmethacrylate  (PMMA) , and  an  analytical 
expression  for  the  flame  spreading  velocity  involving  one  empirical 
constant  was  produced  Agreement  with  experimental  obser  aticn  is 
quite  good,  although  the  authors  caution  the  use  of  this  model  to 
predict  flame  spreading  velocities  tor  situations  in  which  the  fuel 
is  other  than  PMMA. 

Finally,  Feng  and  Sirignano  [25]  have  analytically  treated 
the  process  ol  flame  spreading  over  the  suriace  of  materials  which 
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undergo  an  exothermic  surface  reaction.  Effects  of  fuel  bed  thickness 
and  motion  of  the  gaseous  environment  on  the  flame  spreading  velocity 
were  considered;  conduction  through  the  solid  was  assumed  to  be  the 
dominant  heat  transfer  mode.  However,  since  the  exothermic  reaction 
occurring  in  the  comburtion  of  the  materials  studied  experimentally 
[14-20]  is  not  a surface  reaction  but  rather  takes  place  in  the  gas 
phase,  a meaningful  comparison  of  the  predictions  of  tills  model  with 
experimental  observation  was  precluded. 

In  the  case  of  flame  spread  over  the  surface  of  a liquid 
below  its  flash  point,  the  flame  spreading  velocity  has  been  observed 
to  be  controlled  by  the  flow  of  the  liquid  caused  by  surface  tension 
and  buoyancy  effects.  The  occurrence  of  this  liquid  flow  phenomenon 
has  been  pointed  out  by  the  work  summarized  in  a paper  by  Akita 
[26],  and  by  a more  recent  experimental  study  bv  Helmstetter,  Dryer, 
and  Classman  [27],  Thus,  an  additional  level  of  complexity  is  present 
over  the  case  in  which  the  fuel  is  a solid,  and  attempts  at  analysis 
to  date  have  met  with  only  limited  success  It  is  interesting,  to 
note  that  several  of  the  analytical  efforts  previously  discussed 
on  the  flame  spread  over  solid  surfaces  [21,  23,  24]  originally 
purported  to  be  applicable  to  liquid  surfaces  as  well.  However, 
tue  phenomenon  of  liquid  flow  during  the  flame  spreading  process 
was  not  considered  in  these  models,  and  the  res. .'its  aie  there)  ore 
not  germane  in  this  situation. 

1.2.4  Exothermic  Alloying  of  Metals 

In  the  process  in  which  metals  are  alloyed  to  produce 
enhanced  physical  properties  or  reduced  cost,  n knowledge  of  the 
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time  required  for  the  alloying  reaction  to  take  place  la  of  great 
Importance.  This  fact  has  prompted  several  studies  (a.g,,  28,  29, 
30)  of  the  rates  of  formation  of  various  alloys.  The  particular 
situation  in  which  the  alloying  reaction  is  exothermic  and  self- 
propagating  has  been  Investigated  analytically  by  Hardt  and  Phung 
[31],  and  experimentally  by  Hardt  and  Holslnger  [32].  Self- 
propagating  reactions  were  found  to  take  place  only  in  cases  in 
which  the  reaction  products  were  liquids.  A model  based  on  a mass 
diffusion-controlled,  condensed  phase  exothermic  chemical  reaction 
occurring  ir.  "diffusion  cells"  of  linear  geometry  was  developed, 
and  the  resulting  set  of  equations  was  solved  numerically  in 
dimensionless  form.  A highly  simplified  explicit  solution  was  also 
presented.  Reasonable  agreement  between  analytical  prediction  and 
experimental  measurement  of  reaction  propagation  velocities  was 
obtained. 


1.2.5  Combustion  of  Magnesium-PTFE  Pyrotechnic  Pellets 

In  view  of  a relatively  high  energy  density,  and  the  fact 
that  the  products  of  combustion  are  condensable,  magneslum-PTFE 
pyrotechnic  pellets  have  seen  increasing  use  in  recent  years  as 
ignition  materials  for  use  in  solid  propellant  rocket  engines  and 
similar  applications.  Fabricated  by  compressing  a mixture  of  very 
small  magnesium  particles  and  PTFE  resin,  these  macros cop lcally 
homogeneous  pellets  are  available  commercially  in  a variety  of  sizes 
and  shapes  [33].  In  the  only  work  reported  in  the  literature, 
Griffiths,  O’Sullivan,  and  Thackery  [34]  investigated  flame  shapes, 
flame  spectra,  and  burning  times  during  the  lengthwise  combustion 
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of  cylindrical  magnesium-PTFE  pellets  at  constant  pressure  over  the 
range  0.02-1.0  bar.  Pellets  containing  30  to  60  weight  percent 
magnesium  were  studied,  the  stoichiometric  composition  being  32 
weight  percent  magnesium,  as  given  by  the  overall  combustion  reaction 

Mg  + j C2F4  * MgF2  + C , (1.3) 

which  was  confirmed  by  calorimetry.  The  flame  was  observed  to  be  a 
thin  disc  located  on  the  burning  surface  cf  the  pellet  during 
combustion  at  1.0  bar;  at  lower  pressures,  this  zone  was  seen  to 
expand  and  move  away  from  the  pellet  surface  S-shaped  pressure 
versus  burning  time  curves  were  reported  for  pellets  containing 
50  weight  percent  magnesium,  from  which  the  authors  postulated  the 
existence  of  two  burning  regimes  In  both  cases,  the  burning  rate 
is  assumed  to  be  controlled  by  the  rate  at  which  the  solid  PTFE  is 
thermally  decomposed  to  gaseous  which  then  reacts  with  the 

magnesium.  In  the  high  pressure  burning  regime,  it  is  proposed 
that  the  heat  required  to  thermally  decompose  the  PTFE  comes  from 
the  flame,  while  at  lcwer  pressures,  it  is  reic  that  this  heat  comes 
from  an  exothermic  surface  reaction  No  analytical  mcdel  of  the 
actual  combustion  process  was,  however,  presented 

1 . 3 Specific  Research  Objectives 

In  light  of  the  promising  results  f’-cm  the  preliminary 
experiments  previously  described,  the  chemical  reaction  of  lithium 
with  PTFE  appeared  as  an  interesting  and  i-iatle  high  energy  density 
heet  source  However,  no  previous  work  directly  applicable  to  this 
system.,  either  cf  an  analytical  or  an  experimental  nature,  was 
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reported.  In  order  to  gain  more  insight  into  the  lithium-PTFE 
reaction  process,  then,  further  work  was  indicated-  The  specific 
objectives  of  this  investigation  are,  therefore,  as  follows: 


j 

1 

5 

1.  To  design  and  construct  an  apparatus  with  which  the 

i 

lithium-PTFE  reaction  propagation  process  might  ko  j 

1 

experimentally  examined.  ; 

2.  To  conduct  an  experimental  investigation  of  the  j 

3 

influence  of  varying  the  diameter  and  density  of  PTFE  ■ 

on  the  reaction  propagation  velocity 

3.  To  develop  a model  for  the  lithium-PTFE  reaction 
propagation  process. 

i 

j ;j 
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CHAPTER  II 


THEORETICAL  CONSIDERATIONS 

2 . 1 Description  of  the  Proposed  Reaction  Mechanism 

Thennochemical  calculations  previously  described  have 
indicated  that  the  products  of  the  lithium-PTFE  reaction,  under  the 
conditions  of  interest,  are  in  the  condensed  phase,  Supporting 
evidence  for  this  conclusion  may  be  found  in  the  results  of  the 
preliminary  experiments  described  in  Section  1.1,  in  which  no 
indication  of  high  pressure  within  the  sample  container  tube,  after 
the  reaction  had  taken  place,  was  seen.  Further,  if  it  is  assumed  that 
the  reaction  products  are,  indeed,  gaseous,  and  that  they  are  forced 
to  occupy  the  same  volume  as  originally  taken  up  by  the  reactants, 
as  necessarily  must  be  the  case  in  a sample  container  with  minimal 
free  volume,  a perfect-gas  law  estimation  of  the  pressure  required 
to  achieve  the  dictatea  density  at  reasonable  temperatures  indicates 
that  a pressure  le-el  of  the  order  of  3500  bar  is  necessary.  Since 
the  existence  of  a pressure  of  this  order  cf  magnitude  would  surely 
have  manifested  itself  in  the  experiments  by  deformation  rr  rupture 
of  the  sample  container  tube,  it  is  concluded  that  the  products  of 
the  lithium-PTFE  reaction  are,  in  fact,  in  the  condensed  phase  as 
originally  postulated  The  answer  to  this  impcrtant  question  having 
been  satisfactorily  ascertain.'-.,  the  following  mechanism  for  the 
lithium-PTFE  reaction  process  is  proposed: 

The  overall  reaction  rate-controlling  process  is  assumed 
to  be  the  thermal  decomposition  of  PTFE,  as  set  forth  by  Griffiths, 
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ec  al.,  [34]  for  the  combustion  of  magnesium-PTFE  pyrotechnic 
pellets.  It  is  postulated  that  the  PTFE  decomposes  to  give  the 
gaseous  monomer  C^F^,  found  to  be  the  predominant  decomposition 
product  in  vacuum  pyrolysis  studies  [35*37],  according  to  an  Arrhenius- 
type  first-order  rate  law  obtained  from  vacuum  pyrolysis  data  [37,  38]. 
This  8aa  is  assumed  to  chemically  react  immediately  upon  its 

formation  with  liquid  lithium  to  form  condensed-phase  products 
according  to  the  combustion  reaction  previously  set  forth 

6Li  + C,F.  ♦ 4LiF  + Li_C„  (1.1) 

2 4 2 2 

Heat  transfer  from  the  reacting  lithium  and  forward  to  preheat 

and  thermally  decompose  the  unreacted  PTFE  is  assumed  to  occur  only 
by  condensed-phase  conduction,  and  it  is  by  the  race  at  which  this 
heat  tranafer  takes  place  that  the  reaction  propagation  velocity  is 
controlled . 


2 . 2 One-D. menalonal  Model 

At  the  outset  of  the  analytical  phase  of  this  work,  it 
was  decided  that  a preliminary  investigation  of  the  llchium-PTFE 
reaction  process  at  the  one-dimensional  ie''el  might  prove  instructive. 
It  was  felt  that  even  though  many  of  the  intuitively  important  aspects 
of  this  process  would  be  'gnored  by  necessity  in  a one-dimensional 
description,  an  indication  of  the  validity  of  the  pivotally  important 
assumptions  which  were  made  relative  to  the  identity  of  the  overall 
reaction  rate-controlling  process  and  the  predominant  mcde  cf  heat 
transfer  would  be  gained.  Accordingly,  the  one-dimensional , 
constant  property  form  of  the  heat  conduction  equate cn  with  heat 
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generation  [39], 


|1 . a i_I  + -3_ 

57  3x2  cCP 


(2-1) 


was  solved  numerically  in  dimensional  form,  Stcichicmetric  mixtures 
of  lithium  and  PTFE  only  were  considered,  the  presence  of  lithium 
in  excess  of  the  stoichiometric  amount  and  the  accompanying  loss  or 
heat  to  this  excess  lithium  being  ignored.  The  rate  of  heat  generation 
was  assumed  to  be  determined  by  the  rate  of  thermal  decomposition  of 
the  PTFE,  and  the  conduction  of  heat  through  the  PTFE  was  .onsldered 
negligible.  Although  several  additional  highly  unrealistic 
assumptions  were  necessarily  made  out  of  deference  tc  the  desire  for 
a one-dimensional  representation,  the  predicted  reaction  propagation 
velocity  was  approximately  240  m/s,  which  is  of  the  same  order  cf 
magnitude  as  the  reaction  propagation  vei; cities  measured  in  preliminary 
experiments.  The  fact  that  the  predicted  velocity  was  ct  the  correct 
order  of  magnitude  was  considered  highly  encouraging  It  was  felt 
that  the  two-dimensional  effects  tc  be  m .1. dad  m s_bsequer. : 
models  would  produce  a decrease  in  the  ptediccec  reaction  propagation 
velocity,  bringing  it  into  better  agreement  with  experirentai 
observations  - 


2 3 Twc -Dimensional  Model 

The  two-dimensional  model  cf  the  lithium-FTFE  reaction 
process  is  shown  schematically  in  Figure  3 As  seen  in  this  sketch, 
the  model  is  in  ta'tesian  coordinates , wich  a aye c of  PTFE  or 
thickness  located  adjacent  to  a lays1'  :i  ixthi.T  cf  total  thickness 
The  amount  of  lithium  required  tc  react  s tcithiometrically 

Best  Available  Cor 


Figure  3 Lithium-PTFE  Reaction  Model 
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with  the  PTFE  is  of  thickness  h^g,  the  excess  lithium  layer  is  of 
thickness  h^,  snd  the  length  of  the  sample  is  £.  The  reaction  tone, 
which  moves  from  left  to  right,  contains  unreacted  lithium, 
unreacted  PTFE,  and  reaction  products.  Heat  generated  within  the 
reaction  zone  is  transferred  by  conduction  to  the  lithium,  to  the 
PTFE,  and  to  the  reaction  products.  As  shown  m Figure  3,  the  reaction 
products  are  assumed  to  be  in  the  condensed  phase,  and  to  occupy 
the  volume  originally  taken  up  by  the  PTFE  and  the  stoichiometric 
lithium.  The  llthium-PTFE  reaction  process  is  assumed  to  be  initiat'd 
by  the  action  of  an  ignition  source,  which  serves  to  start  the  chemical 
reaction  of  a small  amount  of  lithium  and  PTFE  located  at  the  x»0 
end  cf  rhe  sample  adjacent  to  the  lithium-PTFE  interface,  and  to 
maintain  the  reaction  products  produced  at  this  location  at  the 
stoichiometric  adiabatic  flame  temperature  It  is  further  assumed 
that  there  are  no  external  heat  losses  from  the  sample. 

Since  the  overall  reaction  rate  is  assumed  to  be  controlled 
by  the  rate  of  thermal  decomposition  of  the  PTFE,  it  is  therefore 
postulated  that  the  reaction  actually  occurs  a;  the  decomposing 
surfa;e  of  the  PTFE.  For  this  region,  the  twc-dimensicnal  constant 
property  form  of  the  heat  conduction  equation  with  he=»t  generation  [39 j. 


may  be  written.  In  the  stoichiometric  lithium  layer  and  in  the 
excess  lithium  layer,  where  nc  reaction  is  assumed  t:  take  place,  the 
temperature  field  may  be  described  by  the  following  equation  (39): 
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which  is  identical  to  Equation  (2.2)  with  the  exception  that  no  heat 
source  term  is  present.  Although  ccnstnnt  property  forms  of  the  heat 
conduction  equation  have  been  written,  it  is  recognized  that  the 
thermal  diffusivity,  a,  in  the  above  two  equations  is  not  constant, 
but  rather  is  a function  of  the  mass  fractions  of  reactants  and 
products  present  in  the  PTFE  and  the  stoichiometric  lithium  liyers. 
This  fact  will  be  accounted  for  in  the  subsequent  solution  of  these 
equations. 

In  order  to  arrive  at  an  expression  tor  the  rate  of  heat 
generation  per  unit  volume,  <j , occurring  in  the  PTFE  layer,  attention 
is  focused  on  the  Arrhenius-type  first-order  rite  equation  for  the 
vacuum  pyrolysis  of  PTFE  (37), 


dm 

<Ft 


-E  /RT 
, 3 

-A  me 
o 


(2.<0 


where: 

m • undecomposed  mass  of  PTFE. 

The  mass  fraction  of  undecomposed  PTFE,  ZT* , is  new  denned  as 
follows: 


Z 


* 

T 


(2.5) 


where: 


m^  ■ original  mass  of  PTFE. 

Substitution  of  Equation  (25)  into  Equation  (24)  and  cancellation  of 


common  terms  yields: 
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-E  /RT 

-A  Z *e  8 
o T ® 


(2.6) 


An  expression  for  the  heat  generation  per  unit  volume  occurring  in 
the  PTFE  layer  may  then  be  written  as  follows: 

dZ  * 

« - -W  -jr  • <2  ’> 


the  negative  sign  arising  since  the  rate  ot  heat  gene  rati :n  is 
proportionaly  to  the  rate  of  evolution  ot  'apor,  which  is 

numerically  equal  to  the  rate  of  decrease  of  the  mass  of  undecomposed 
PTFE,  but  cf  opposite  sign. 

The  two-dimensional  analytical  model  for  the  lithium-PTFE 
reaction  process  therefore  consists  ot  the  following  set  ot  equations: 
In  the  PTFE: 


where : 


3T 

dt 


dZ?* 

dt 


/W 

|dY 

(2.8) 

1 dx2  dy2/ 

1 ^CP 

1 dt 

-E  -RT 

■ -A  Z *e  J 
o 1 

(2.6) 

i stoichiometric 

and  exces 

s aichium: 

dT  / 32T  a d2T 

St  ’ aL  2 ,2 

| ux  dy 


(2.3) 


The  appropriate  boundary  and  initial  conditions  are  as  fellows: 
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Boundary  Conditions: 


3T 


x ■ 0,  y f hT  ; - 0 


x - 0,  y - hT  ; T - , 


x - A ; 


£-° 


y-o  , |I-o 
»-■>  > 


Initial  Conditions: 


(2.9) 


t 


t 


t 


0 ; T • To  , all  x,  y,  except 

0 ; T “ Taf  * x “ °*  y “ hT  ' 
0 ; ZT*  • 1 . 


(2.10) 


In  order  to  make  the  model  more  general,  and  to  reduce  the 
number  of  parameters  which  characterize  the  system,  the  following 
dimensionless  variables  are  defined: 


T* 


» 


(2.11) 
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where : 

■ the  stoichiometric  adiabatic  flame  temperature  at  one 
atmosphere  total  pressure 

and 

t ^ ” the  time  required  to  thermally  decompose  an  arbitrary 
fraction  of  an  amount  of  PTFE  at  thiB  stoichiometric 
adiabatic  flame  temperature. 

Substituting  these  dimensionless  variables  into  Equations 
(2.3),  (2.7),  and  (2.8)  yields: 

In  the  PTFE: 


3T*  a [32T*  32T*  l PTAHRW|  dZT* 

3Y*2J"  pCPTafl  dl* 


dZT* 


■<Vch> 


Z *e 


-(Ea/RTaf)  * 


v.2  . XL  ) 


(2.13) 


In  the  stoichiometric  and  excess  lithium: 


3T*  _ o_  32T*  32T* 

3t*  °T  V3X*2  3Y*2 


(2. la) 


The  following  dimensionless  parameters  are  now  defined: 


A*  - — , 

aT 


R* 


ot^hrw 


F*  - A t , 

o ch 


E*  * ITT 


af 


(2.15) 
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The  parameters  F*  and  E*  above  are  constants.  The  parameter  A*  is 
In  effect  a variable,  its  value  depending,  first,  on  whether  it  it 
evaluated  in  the  stoichiometric  lithium  layer  or  in  the  PTFE  layer, 


and  second,  on  the  mass  fractions  of  reactants  and  products  present 
at  the  location  in  equation.  Similarly,  the  value  of  the  parameter 


R*,  which  pertains  only  to  the  PTFE  layer,  is  a function  of  the 
chemical  composition. 

The  dimensionless  analytical  model  is  then  given  by  the 
following  set  of  equations  with  appropriate  boundary  and  initial 
conditions : 


In  the  PTFE: 


. AJihi  + £zl\  . R*  fil 

^ lax*2  3Y*2 1 d’**  ' 


(2.16) 


dZ  * 

ZI_  - _F*  z * -e*/t* 

dr*  * i e 


(2.17) 


In  the  stoichiometric  and  excess  lithium: 


m . AJ  ii ii  + i!n 

\ 3X*  dY* 
Boundary  Conditions: 


(2.18) 


X*  - 0 , Y*  * Hj*  ; 


- 0 , 


X*  - 0 , Y*  - Ht*  ; T*  - 1 


X*  = L* 


’ 3X* 


Y*  = 0 


> riv* 


’ 3Y* 


Y*  - H* 


(2.19) 
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Initial  Conditions: 


T*  • 0 ; 

X*  ■ T * 

0 

, all  X*,  Y*.  except 

T*  • 0 ; 

T*  « 1 

, X*  - 0,  Y*  - , 

- 0 ; 

v-1 

• 

(2.20) 

2,4  Details  of 

the  Numerical  Solution 

As  a first  step  in 
dimensional  model,  attention 

the  numerical  solution  of  the  two- 
was  focused  cpon  the  dimensional 

form 

of  this  model,  as  given  by  Equations  (2.3),  (2.6),  (2  8),  (2.9)  and  (2.10). 
The  partial  differential  equations  were  cast  into  finite  difference 
form,  and  a digital  computer  program  for  the  numerical  solution  of 
the  model  was  written.  The  purpose  of  this  effort  was  to  ascertain 
whether  the  predicted  reaction  propagation  velocities  were  of 
reasonable  magnitude,  and  to  identify  the  variables  to  which  the 
reaction  propagation  velocity  was  most  sensitive.  The  particular 
techniques  used  in  this  numerical  solution  were  identical  to  those 
used  in  the  solution  of  the  dimensionless  form  :r  the  two-dimensional 
model,  which  will  be  described  subsequently  The  validity  of  the 
numerical  solution  cf  the  partial  difterential  equations  was  checked 
by  comparison,  under  appropriate  boundary  conditions  and  with  no 
heat  generation,  with  the  analytical  solution  for  the  fansient 
heating  of  a semi-infinite  slab  [39],  and  was  found  to  agree  to  within 
0.5  percent  The  numerical  solution  of  the  ordinary  differential 
equation  under  the  condition  of  constant  temperature  was  compared  with 
its  analytical  solution  with  similar  results  For  conditions  similar 
to  those  present  in  preliminary  experiments,  the  predicted  reaction 
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propagation  velocities  were  approximately  15X  higher  than  those 
observed,  which  was  considered  to  be  reasonable  agreement-  Further, 
through  a systematic  variation  of  parameters,  it  was  discovered 
that  the  reaction  propagation  velocity  was  most  sensitive  to  small 
changes  in  the  thermal  dlffusivity  of  the  PTFE,  a , and  in  the 
activation  energy  for  the  thermal  decomposition  of  the  PTFE,  E . 

These  results,  then,  provided  the  clue  for  the  characteristic  time 
and  length  scales  of  the  system,  which  were  used  as  normalizing 
parameters  in  the  nondimensionalization  of  the  model  as  given  by 
Equation  (2.11). 

The  dimensionless  form  of  the  two-dimensional  model,  as 
given  hy  Equations  (2.16-2.20),  was  then  considered.  The  llthium-PTFE 
reaction  sample  was  first  broken  up  into  a grid  network  as  shown  in 
Figure  4.  As  seen  in  Figure  4,  the  X*  direction  grid  spacing,  AX*, 
was  the  same  for  the  PTFE  layer  as  for  the  stoichiometric  and  excess 
lithium  layers.  For  the  two  lithium  layers,  the  Y*  direction 
grid  spacing,  AYL*.  was  also  the  same.  The  Y*  direction  grid  spacing 
in  the  PTFE  layer,  AY^,*,  was  however  adjusted  such  thac  the  amount 
of  PTFE  represented  by  each  grid  point  in  this  layer  wa6  just  that 
amount  necessary  to  react  stoichiomecrically  with  the  amount  of 
lithium  represented  by  each  grid  point  in  the  stoichiometric  lithium 
layer.  The  Y*  direction  grid  spacing  spanning  the  lithium-PTFE 
interface,  AY^*,  was  therefore  equal  to  the  average  of  AY^*  and  AY^*. 
The  number  of  grid  points  in  the  PTFE  layer  was  thus  equal  to  the 
number  of  grid  points  in  the  stoichiometric  lithium  layer;  the  number 
of  grid  points  in  the  excess  lithium  layer  was  arbitrary,  but,  once 
chosen,  determined  the  amount  of  excess  lithium  present.  The  two  heat 


Figure  4 Sketch  of  Finite  Difference  Grid  Network 


conduction  equations,  Equations  (2.16)  and  (2.18)  were  cast  in  finite 
difference  form  uoing  the  central  difference  approximation  for  the 
spatial  derivatives  and  the  forward  difference  approximation  for  the 
time  derivatives  (co)  ; the  detailed  nodal  equations  are  presented  In 
Appendix  A.  The  ordinary  differential  equation.  Equation  (2.17) 
was  put  into  form  tor  numerical  Integration  through  the  use  of 
Euler's  method  (41). 

For  the  purposes  of  the  numerical  solution  of  the  model,  it 
was  assumed  that  the  chemical  reaction  between  the  PTFE  and  the 
stoichiometric  lithium  took  place  according  tc  the  following  sequence 
As  seen  in  Figure  4,  the  initial  conditions  given  by  Equation  (2.20) 
dictated  that  reaction  products  at  7*»1  were  presen:  at  grid  points 
1,  nT,  and  1,  nT+l.  The  reaction  zone  was  allowed  to  propagate 
along  the  lithium-PTFE  interface  at  a rate  determined  by  heat 
transfer  and  PTFE  decomposition  kinetics,  the  FTFE  represented  by 
grid  points  i,  nT  reacting  stoichiometncally  with  the  lithium 
represented  by  grid  points  i,  n^+1.  As  the  reaction  zcne  spread 
into  the  PTFE  and  the  stoichiometric  lithium  in  the  direction 
perpendicular  to  the  lithium-PTFE  inte-face,  it  was  assumed  that  the 
PTFE  represented  by  grid  points  i,  n^-1  rea'ted  with  the  .lithium 
represented  by  grid  points  i,  nT+2  This  lateral  spreading  of  the 
reaction  zone  was  allowed  to  continue  in  this  symmetrical  manner 
through  as  many  grid  points  as  were  located  In  the  PTFE  and 
stoichiometric  lithium  until  the  PTFE  represented  by  grid  points 
i,  1 had  reacted  with  the  lithium  represented  by  grid  points  i, 
n'.  At  any  given  X*  location  (any  value  of  : between  1 and  n ), 

^ A 

the  reaction  was  required  to  have  gone  cc  completion  at  the  pair  of 
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reacting  grid  points  closest  to  the  lithium-PTFE  in. erf ace  before  it 
was  allowed  to  progress  to  the  subsequent  reacting  grid  point  pair. 

For  the  purpose  of  property  computations,  it  was  necessary 
to  monitor  the  mass  fractions  of  unreacted  (undecoaposed)  PTFE  and 
reaction  products,  and  unreacted  lithium  and  reaction  products 
present  in  the  PTFE  and  the  stoichiometric  lithium  layers, 
respectively,  as  a function  of  dimensionless  time.  Initially,  the 
mass  fraction  of  unreacted  PTFE,  Z,p*,  at  each  of  the  grid  points 
in  the  PTFE  layer,  with  the  exception  of  grid  point  1,  nT>  was 
unity  and  the  reaction  product  mass  fraction,  Zp*,  was  zero. 

Similarly,  the  initial  mass  fraction  of  unreacted  lithium,  Z^*,  was 
unity  and  the  reaction  product  mass  fraction,  Zp* , was  zero  at  each 
of  the  grid  points  in  the  stoichiometric  lithium  layer  except  at  the 
grid  point  1,  nT+l.  As  the  PTFE  was  decomposed  and  reacted  at  each 
of  the  grid  points  in  the  PTFE  layer,  the  value  of  Z^*  at  each  point 
decreased  according  to  Equation  (2.17).  The  value  of  Z^*  at  each 
corresponding  grid  point  in  the  stoichiometric  lithium  layer  was 
therefore  decreased  as  dictated  by  the  assumed  symmetrical  nature  of 
the  lateral  spread  of  the  reaction  zone.  Frcm  this  symmetry,  and 
from  the  stoichiometric  sizing  of  the  grid  network  previously 
described,  it  is  seen  from  Figure  4 that: 

ZL*(i,  2^+1  - j)  - ZT*(i,  j>  (2.21) 

and 


Zp*U,  2nT+l  - j)  - Zp*(i,  J)  - 1 - ZT*(i,  J)  . 


(2.22) 
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In  order  to  accomplish  ths  numerical  solution  of  the 
dimensionless  form  of  the  two-dimensional  model,  it  was  necessary  to 
compute  the  value  of  A*,  the  dimensionless  thermal  diffusivity  ratio 
appearing  in  Equations  (2.16)  and  (2.18).  A knowledge  of  the  thermal 
diffusivity,  defined  as 


PC, 


(2.23) 


of  the  lithium,  the  PTFE,  and  the  reaction  products  was  therefore 
required.  The  values  of  the  thermal  conductivity,  density,  and 
specific  heat  used  to  determine  the  thermal  dif f usivities  of  the 
lithium  and  the  PTFE  were  obtained  from  Reference  [42),  and 
Reference  [43],  respectively . The  reaction  products  were  assumed 
to  possess  the  thermal  conductivity  [44],  and  the  specific  heat  [1] 
of  liquid  lithium  fluoride,  which  accounts  for  approximately  73 
weight  percent  of  the  products  according  to  the  postulated  reaction 
stoichiometry.  The  density  of  the  reaction  products  was  determined 
by  the  earlier  assumption  that  these  products  are  in  the  condensed 
phase,  and  fill  the  volume  originally  occupied  by  the  unreacted 
lithium  and  PTFE. 


The  thermal  diffusivity  of  all  substances  is  a function  of 
temperature,  the  value  of  a of  pure  lithium,  tci  example,  increasing 
by  approximately  a factor  of  2 from  the  melting  point  to  the 
normal  boiling  point  142].  However,  during  the  course  cf  the 
postulated  reaction,  in  which  the  stoichiometric  lithium  and  the 
PTFE  are  transformed  into  reaction  products,  the  thermal  diffusivity 
of  the  stoichiometric  lithium  is  seen  to  decrease  by  a factor  of  44, 
and  that  of  the  PTFE  to  increase  by  a factor  of  approximately  6. 
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For  this  re *8 on,  average  values  of  the  thermal  conductivities, 
densities,  end  specific  heats  of  the  lithium,  PTFE,  and  reaction 
products  were  used  in  the  computation  of  the  desired  thermal  dlffuslvlty 
ratios,  temperature  variations  in  these  ratios  being  neglected  with 
respect  to  variations  caused  by  changes  in  chemical  composition. 

In  order  to  compute  the  value  of  the  dimensionless  parameter 
A*  in  the  stoichiometric  lithium  layer  as  a function  of  chemical 
composition,  the  lithium  and  reaction  products  were  assumed  to  form 
a perfect  solution  whose  thermal  conductivity,  specific  volume,  and 
specific  heat  were  linear  functions  of  mass  fraction.  The  thermal 
dlffuslvlty  ratio  A*  for  this  layer  is,  therefore,  given  by  the 
following  relationship: 


A* 


•s  kp  1/  P-j.  Pf  1 
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(2.24) 


The  PTFE,  however,  was  assumed  to  retain  Its  original  properties 
until  completely  reacted.  As  the  reaction  progressed  at  each  grid 
point  In  the  stoichiometric  lithium  layer,  then,  the  value  of  A*  was 
computed  from  Equation  (2.24),  At  the  corresponding  grid  point  in 
the  PTFE  layer,  the  value  of  A*  was  held  et  unity  until  the  reaction 
had  gone  to  completion,  at  which  time  A*  was  set  equal  to  that  for  the 
reaction  products: 


A*  ■ — — 

aT 


(2,25) 
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Consistent  with  the  assumption  that  the  PTFE  retained  its 
original  properties  until  reacted,  the  values  for  the  density  and 
specific  heat  of  PTFE  were  used  in  the  denominator  of  the  dimensionless 
parameter  R*,  as  given  by  Equation  (2.15), 

In  order  to  simplify  the  numerical  eor-icion  ::  the  model, 
the  melting  processes  occurring  in  the  lithium  and  in  the  PTFE,  and 
the  heat  absorbed  in  the  vaporization  or  the  PTFE  were  ignored.  It 
was  anticipated  that  erroneously  high  temperatures  might  be  predicted 
by  the  numerical  solution  of  the  model  as  a result  of  neglecting 
these  heat-absorbing  phenomena,  as  well  as  from  the  use  at  average 
heat  transfer  properties  as  previously  described,  The  predicted  value 
of  T*  in  the  numerical  solution  was  therefore  limited  to  unity, 
corresponding  to  the  assumption  that  the  temperature  at  any  point 
within  the  reaction  sample  could  not  exceed  the  stoichiometric 
adiabatic  flame  temperature  for  one  atmosphere  total  pressure.  Under 
actual  experimental  conditions,  temperatures  greater  than  the 
stoichiometric  siiabifc  flame  tempe1- atu-e  a*  standa-d  conditions 
are  possible  due  to  the  energy  addition  tc  the  reaction  sample 
caused  by  the  firing  of  the  igniter,  and  the  la-k  :i  assurance  that 
the  reaction  dees  not  take  place  at  presences  in  excees  of  one 
atmosphere.  However,  in  reality  there  9re  also  heat  losses  from  the 
sample,  which  ha”e  an  opposing  effect  In  light  of  these  facts, 
then,  it  is  felt  that  this  assumption  is  not  unreasonable. 

The  nature  of  the  solution  or  the  PTFE  chermal  decomposition 
rate  law,  gi’'en  in  dimensionless  form  by  Equation  (2  i’\  is 
logarithmic,  the  value  of  Z^*  approaching  zero  aeympt -t: tally . Since 
an  infinite  amount  of  dimensionless  time  would  be  required  tor  the 
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undecomposed  PTFE  mass  fraction  to  exactly  go  to  tero,  the  chemical 
reaction  at  grid  points  in  the  PTFE  layer  was  assumed  to  have  gone 
to  completion  when  the  value  of  dropped  below  the  value  0.10, 
corresponding  tc  90%  of  the  PTFE  having  been  decomposed  and  chemically 
reacted. 

A digital  computer  program  to  accomplish  the  numerical 
solution  of  the  dimensionless  form  of  the  two-dimensional  model  in 
the  manner  described  above  was  written  in  Fortran  IV  language.  A flow 
diagram  of  the  computer  program  is  presented  in  Appendix  B-  The 
numerical  solution  of  the  dimensionless  heat  conduction  equations, 
Equations  (2.16)  and  (2.18),  with  no  heat  generation  and  under 
appropriate  boundary  conditions,  was  again  compared  with  the 
analytical  solution  for  the  transient  heating  of  a semi-infinite 
slab  [39],  and  was  again  found  to  agree  to  within  0.3  percent.  The 
numerical  solution  of  the  dimensionless  PTFE  thermal  decomposition 
rate  equation.  Equation  (2.17),  was  also  checked  under  constant 
temperature  conditions  with  its  closed-fora  solution  with  similar 
results.  In  addition,  an  energy  balance  check,  in  which  the  energy 
absorbed  in  the  heating  of  the  reactants  and  products  was  compared 
to  the  energy  liberated  by  the  chemical  reaction,  was  incorporated 
into  the  computer  program  to  insure  that  the  numerical  simulation 
was  thermodynamically  correct 

The  dimensionless  reaction  propagation  velocity  was  determined 
by  continuously  monitoring  the  product  mass  fraction  at  the  grid 
points  in  the  stoichiometric  lithium  and  in  the  PTFE  immediately 
adjacent  to  the  lithium-PTFE  interface  As  the  reaction  zone  moved 
along  the  lithium-PTFE  interface  and  reaction  was  completed  at  each 
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successive  pair  of  grid  points  adjacent  to  the  interface,  the 
dimensionless  reaction  propagation  velocity  was  computed  by  dividing 
the  dimensionless  distance  traversed,  AX*,  by  the  dimensionless 
time  elapsed  since  the  previous  pair  of  grid  points  had  completely 
reacted.  Examination  of  the  dimensionless  X*  and  t*  variables  as 
given  by  Equation  (2.11)  then  yields  the  following  definition  of  the 
dimensionless  reaction  propagation  velocity: 


V* 


(2.26) 


where 


2 . 5 Results  of  the  Analytical  Model 

For  the  numerical  solution  of  the  dimensionless  form  of  the 
two-dimensional  model,  the  values  of  the  thermal  conductivity, 
density,  and  specific  heat  of  the  lithium  were  taken  from  Reference 
[42]  as  those  for  liquid  lithium  at  800  K.  The  properties  of  the  FIFE 
were  obtained  from  Reference  [43],  while  the  values  ox  the  frequency 
factor  and  the  activation  energy  for  the  thermal  decomposition  of 
PTFE,  as  used  in  Equation  (2.13),  were  taken  from  Reference  [37],  The 
thermal  conductivity  and  specific  heat  of  the  reaction  products 
were  obtained  from  Reference  144]  and  Reference  [1],  respectively, 
as  chose  for  liquid  lithium  fluoride  at  1500  K,  the  density  of  the 
reaction  products  being  determined  as  previously  described.  The 
values  of  the  physical  and  chemical  properties  used  in  the  numerical 
solution  of  the  model  are  given  in  Appendix  C 
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Using  molecular  weight  dace  from  Reference  [1],  the 
stoichiometric  PTFE-lithium  mass  ratio  was  calculated  from  the  assumed 
overall  reaction  mechanism,  as  given  by  Equation  (1.1).  The  standard 
heat  of  formation  of  solid  PTFE  was  determined  by  combining  the 
standard  heat  of  formation  of  gaseous  C^F^  obtained  from  Reference  [1] 
with  the  enthalpy  of  depolymerisation  of  PTFE  as  given  in  Reference 
[37).  This  quantity  was  then  used  in  conjunction  with  additional 
data  from  Reference  [1]  to  calculate  the  standard  heat  of  reaction 
of  lithium  with  PTFE,  and  the  stoichiometric  adiabatic  flame 
temperature  of  this  reaction  at  one  atmosphere  pressure.  The 
stoichiometric  adiabatic  flame  temperature  under  standard  conditions 
was  computed  to  be  1966  K,  the  normal  boiling  point  of  lithium 
fluoride,  and  the  value  of  the  characteristic  time,  t ^ , used  as  a 
normalizing  parameter  in  Equation  ''2.11)  was  taken  to  be  the  time 
required  to  thermally  decompose  901  of  an  initial  amount  of  PTFE  at 
this  temperature. 

The  propagation  of  the  reaction  zone  through  the  lithium- 
PTFE  reaction  sample,  as  predicted  by  the  numerical  solution  of  the 
model,  occurs  in  the  following  manner:  Ac  the  start  of  the  solution, 

a period  of  time  elapses  during  which  n:  reaction  occurs.  In  this 
"induction  period,"  heat  is  conducted  from  the  ignition  source, 
located  at  the  X*»0  end  of  the  sample  at  the  lithium-PTFE  interface, 
into  the  lithium  and  PTFE,  preheating  the  reactants.  Following 
this  period,  the  reaction  zone  starts  to  move  through  the  sample, 
experiencing  an  initial  rapid  acceleration.  As  the  reaction  zone 
approaches  and  passes  through  the  center  of  the  sample,  the 
dimensionless  reaction  propagation  velocity  approaches  a constant 
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value,  this  value  rapidly  increasing  once  again  as  the  reaction 
zone  nears  the  X*«l*  end  of  the  sample.  For  values  of  the  dimensionless 
parameters  A*,  R*,  F*,  and  E*  determined  as  previously  described,  and 
a total  PTFE  to  lithium  mass  ratio  of  0.379,  the  predicted  motion  of 
the  reaction  zone  is  illustrated  in  Figure  5-  In  Figure  5,  the 
dimensionless  reaction  propagation  velocity  V*  is  plotted  as  a 
function  of  dimensionless  distance  X*  for  two  values  of  the 
dimensionless  sample  aspect  ratio,  L*/H*-  As  seen  from  this  plot, 
the  solutions  are  essentially  Identical  up  to  the  point  X**3.0. 

Beyond  this  point,  the  dimensionless  reaction  propagation  velocity 
becomes  move  and  more  a function  of  L*/H*,  indicating  that  the 
propagation  of  the  reaction  zone  is  strongly  effected  by  the  presence 
of  the  ends  of  the  sample.  At  L*/H*“3,  these  end  effects  are  quite 
pronounced,  the  value  of  V*  increasing  over  the  entire  length  of  thfc 
sample.  At  L*'H*»10,  however,  the  dimensionless  reaction  propagation 
velocity  is  seen  tc  reach  a level  which  is  essentially  constant 
ever  the  central  portion  of  the  sample  Were  it  possible  to  increase 
the  value  of  the  dimensionless  sample  aspect  ratio  substantially, 
it  is  believed  that  the  predicted  dimensionless  reaction  propagation 
velocity  would  indeed  reach  a rigorously  constant  level.  Computer 
storage  capacity  and  run-time  considerations,  however,  made  this 
prohibitive . 

One  very  important  result  of  this  analytical  investigation 
was  the  predicticn  of  the  steady  dimensionless  reaction  propagation 
velocity  as  a function  of  the  tccai  PTFE-ii'.hium  mass  ratio.  Based 
on  the  assumed  overall  combustion  reaction  given  in  Equation  (1.1), 


Figure 
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Che  stoichiometric  PTFE-lithium  mass  ratio  is  2.402-  In  order  to 
limit  the  final  equilibrium  temperature  of  the  material  in  the  reaction 
sample,  after  the  reaction  had  taken  place,  to  a level  below  the 
melting  point  of  standard  materials  of  construction  such  as  stainless 
steel,  lithium  in  excess  of  that  required  to  react  stoichiometrically 
with  the  PTFE  must  be  present.  Eased  on  the  standard  heat  of  formation 
of  solid  PTFE  determined  as  previously  described,  and  data  from 
Reference  [1],  the  dependence  of  the  final  equilibrium  temperature 
on  the  PTFE-lithium  mass  ratio  is  shown  in  Figure  6.  From  Figure  6 
it  is  seen  that,  in  order  to  avoid  excessive  final  equilibrium 
temperatures,  a maximum  PTFE-lithium  mass  ratio  oi  approximately 
0.4  is  indicated,  the  discontinuity  in  this  plot  being  due  to  the 
heat  absorbed  during  the  fusion  of  lithium  fluoride  at  1121  K. 

The  lower  limit  of  the  PTFE  lithium  mass  ratio  range  to  be 
analytically  investigated,  as  dictated  by  computer  storage  and  run-time 
considerations,  was  0,2,  corresponding  to  a final  equilibrium 
temperature  of  980  K 

As  seen  in  Figure  5,  the  dimensionless  reaction  propagation 
velocity  does  not,  in  fact,  read;  a truly  steady  level  for  either  of 
the  dimensionless  sample  aspect  ratios  examined,  although  the  numerical 
solution  for  L*/H*«10  predicts  a velocity  which  is  very  nearly  steady. 
It  was  thus  decided  to  investigate  the  variation  in  the  dimensionless 
reaction  propagation  velocity,  observed  tc  be  steady  to  within  some 
tolerance  level,  as  a function  of  total  PTFE-lithium  mass  ratio  for 
a dimensionless  sample  aspect  ratio,  L*/K*,  of  10  The  criterion 
for  the  selection  of  the  steady  value  of  V*  was  taken  to  be  the 
point  at  which  the  value  of  v*  was  constant  with  X*  to  within  1%.  A 


UAUUL  .''VIA-  ifal  Aiial  tl 


ploc  of  che  dimensionless  steady  reaction  propagation  velocity 
determined  in  this  manner,  Vg*,  as  a function  of  the  total  PTFE- 
llthium  mass  ratio  is  shown  in  Figure  7.  As  seen  from  this  plot, 
the  value  of  Vg*  decreases  linearly  with  increasing  PTFE-lithium 
mass  ratio  over  the  range  Investigated.  Over  this  range,  the 
relationship  between  Vg*  and  the  PTFE-lithium  mass  ratio  is: 

Vg*  ■»  1,8575  - 1.1536  M*  , (2.28) 

where 


M*  ■ PTFE-lithium  mass  ratio. 

In  order  to  illustrate  the  predicted  dimensionless  temperature 
profiles  present  in  the  reaction  sample  during  the  propagation  of 
the  reaction  front,  It  was  necessary  to  obtain  a numerical  solution  of 
the  Kodei  for  a reaction  sample  with  a small  aspect  ratio,  L*/H*. 

This  adjustment  of  the  aspect  ratio  was  required  in  order  that  the 
desired  information  could  be  represented  with  reasonable  detail  on  a 
sketch  of  manageable  physical  size.  The  predicted  dimensionless 
temperature  profiles  for  a reaction  sample  with  a value  of  L*/H*  of 
1.7,  and  a PIFE-lithiuir  mass  ratio  of  0,206,  are  shown  in  Figure  S 
In  this  sketch,  the  reaction  front  is  seen  to  have  moved  approximately 
1/3  of  the  distance  through  the  sample,  the  dimensionless  reaction 
propagation  veloc.  try  at  tins  point  being  constant  with  X*  to  within  2.5%. 
L\  Figure  8,  ’hie  i sc  them  present  at  the  ignition  source  location 

may  be  seen  in  addition  to  the  T*«l  isotherm  lO'  sted  in  the  reaction 
zone.  Two  f**0  9 isotherms,  one  resulting  from  the  influence  of  the 
ignition  source,  and  the  other  from  the  influence  c l the  heat 
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generated  within  the  reaztion  zone,  are  also  present.  The  temperature 
gradients  within  the  PTFE  are  aeen  to  be  very  much  steeper  than  those 
within  the  lithium  and  reaction  products.  In  addition,  considerable 
heat  conduction  from  the  reaction  products  through  the  excess  lithium 
to  the  PTFE  in  the  region  ahead  of  the  reaction  zone  is  evidenced. 

For  the  numerical  solutions  of  the  dimensionless  form  of  the 
two-dimensional  model  described  above,  the  value  of  the  unreacted 
PTFE  mass  fraction  Z^*,  at  which  the  reaction  was  assumed  to  have 
gone  to  completion  was  taken  as  0.1,  corresponding  to  a product 
mass  fraction  of  0.9.  Although  this  implies  that  only  902  of  the 
original  mass  of  PTFE  present  at  each  grid  point  in  this  layer  was 
taken  into  account  in  the  solutions,  a decrease  in  this  cutoff 
value  to  0.05,  corresponding  to  952  of  the  PTFE  mass  reacted,  had  the 
effect  of  increasing  the  predicted  12  tolerance  steady  value  of  the 
dimensionless  reaction  propagation  velocity  by  less  than  one  percent. 

A PTFE  mass  fraction  cutoff  value  of  01  was,  therefore,  used  in 
these,  as  well  as  in  subsequent,  numerical  solutions. 

As  previously  indicated,  the  dimensionless  form  of  the 
PTFE  thermal  decomposition  rate  law.  Equation  (2.17),  was  solved 
through  the  use  of  Euler's  method  [41],  Since  this  is  a somewhat 
crude  numerical  representation  of  the  solution  cf  an  ordinary 
differential  equation,  the  method  of  Runge  and  Kutta  [41]  was  also 
tried  for  comparison.  Although  the  Runge-Kutta  method  is  considerably 
more  sophisticated,  involving  four  additional  calculations  at  each 
increment  of  che  independent  variable  and  a correspondingly  longer 


total  computation  time,  the  resulting  solution  was  seen  to  be 
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Indistinguishable,  to  five  significant  figures,  from  that  produced 
by  Euler's  method  for  independent  variable  increment  sizes  of 
interest.  Since  it  was  desired  that  computation  time  be  minimized, 
Euler's  method  was,  therefore,  used  in  the  numerical  solution  of  the 
dimensionless  form  of  the  two-dimensional  model. 


CHAPTER  III 


EXPERIMENTAL  INVESTIGATIONS 

3 . 1 Experimental  Apparatus 

The  test  apparatus  which  contained  the  lichlum-PTFE 
reaction  sample  is  shown  in  Figure  9.  Thi6  device  was  similar  in 
most  respects  to  the  improved  test  apparatus  used  in  the  preliminary 
experiments  previously  described.  The  major  component  of  the  test 
apparatus  was  the  sample  tube,  consisting  of  a 216  mm  long,  28.6 
mm  outside  diameter  304  stainless  steel  tube  having  0.89  mm  thick 
walls.  The  igniter  plate,  which  was  furnished  with  a tapped  hole 
to  accommodate  the  igniter,  was  fusion  welded  into  one  end  of  the 
sample  tube.  This  component  was  constructed  of  304  stainless  steel 
as  were  the  remaining  elements  of  the  test  apparatus.  The  cap 
incorporated  redundant  high  temperature  Viton  0-rings  to  form  a 
seal  at  the  end  of  the  sample  tube  opposite  the  igniter  plate,  The 
1 6 mm  thick  wafer  was  placed  in  the  end  of  the  sample  tube  adjacent 
to  the  cap  in  order  to  inhibit  the  flow  of  hot  lithium  and  reaction 
products  intc  the  cap  and  associated  plumbing  after  the  reaction  had 
taken  place.  A small  quantity  of  lithium  fluoride  was  also  placed 
next  to  the  wafer  at  the  end  of  the  PTFE  sample  in  order  to  assist 
in  this  effort  The  wafer  was  provided  with  several  1.6  mm  diameter 
through-holes  to  allow  pressure  communication  with  the  interior  of 
the  sample  tube  The  retainer,  located  at  the  igniter  end  of  the 
sample  tube,  was  so  constructed  as  to  allow  access  to  the  igniter 
(Holex  Model  6102)  while  the  retainer  was  in  place.  Four  tie  rods, 
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two  of  which  ara  shown  in  Figure  9,  were  used  in  conjunction  with  the 
retainer  to  hold  the  cap  in  placa  over  the  end  of  the  sample  tube. 

The  chromel-alumel  thermocouple  probes,  having  1,0  mm  diameter 
sheaths,  0.15  mm  diameter  thermocouple  wire,  and  exposed  measuring 
Junctions,  were  constructed  from  "Omegaclad"  metal  sheathed  thermocouple 
wire.  The  probes  were  installed  through  the  wall  of  the  sample  tube 
at  51  mm  Intervals,  the  thermocouple  measuring  junctions  being  located 
at  the  lithium-PTFE  Interface.  The  thermocouple  probes  were  held  in 
place  by  high  temperature  epoxy  furnished  by  Omega  Engineering, 
Incorporated . 

In  order  to  preclude  problems  should  a high  temperature 
lithium  leak  develop,  the  test  apparatus  was  enclosed  in  a safety 
vessel  while  testing  was  ir.  progress.  The  apparatus  was  mounted 
to  the  removable  lid  of  the  safety  vessel,  which  was  provided  with 
feed-throughs  for  electric  power  and  instrumentation  signals  along 
with  tube  fittings  to  allow  communication  with  the  sample  tube  and 
the  interior  of  the  safety  vessel.  A photograph  of  the  test  apparatus 
attached  to  the  lid  of  Che  safety  vessel  is  shown  cr.  Figure  10.  The 
thermocouple  leads  are  seen  to  have  been  connected  to  extension 
wires  which  terminate  at  a thermocouple  jack  panel  also  mounted  to 
the  safety  vessel  lid.  The  igniter  power  lead  connected  to  the 
igniter  at  the  top  of  the  test  apparatus  is  also  visible  in  Figure  10. 

A sketch  of  the  overall  test  facility  is  shown  in  Figure 
11.  A solenoid  valve  was  included  in  the  test  apparatus  plumbing 
in  order  to  allow  isolation  of  the  reaction  sample  prior  to  the 
initiation  of  each  test  run.  Through  the  use  of  the  argon  tank,  tne 
vacuum  pump,  and  the  3-way  valves  and  associated  plumbing,  it  was 
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Figure  11  Sketch  of  lest  Facility 


56 


possible  co  evacuate  or  pressurize  Che  reaction  sample  and  the  interior 
of  the  safety  vessel  as  dasired,  pressure  gauges  being  provided  to 
monitor  the  pressure  levels  present  at  these  locations. 

A photograph  of  the  test  facility  is  shewn  in  Figure  12. 

The  safety  vessel  and  its  lid  with  the  test  apparatus  attached  are 
visible  at  the  left  of  this  photograph.  The  ccntrcl  panel,  on  which 
were  mounted  the  pressure  gauges  and  associated  plumbing,  and  the 
electric  circuitry  which  controlled  the  operation  of  the  solenoid 
valve  and  che  igniter,  may  also  be  seen  at  the  right  of  the  photograph. 
Time-referenced  records  of  the  temperature  levels  present  at  the 
measuring  junctions  of  the  thermocouple  probes  during  test  runs  were 
produced  by  recording  the  thermocouple  analog  voltage  output  signals 
on  a Tektronix  Model  564  dual  trace  storage  oscilloscope  and  on  a 
CEC  Model  5-124  recording  oscillograph.  The  voltage  signals  from  all 
three  thermocouples  located  at  the  lithiun-PTFE  interface  in  the 
reaction  sample  were  recorded  by  the  oscillograph,  while  only  the 
thermocouple  nearest  the  igniter  and  the  thermo :ouple  nearest  the 
cap  were  connected  to  the  cscll lose epe ; the  oscillograph  and  the 
oscilloscope  are  visible  on  the  table  in  the  center  or  Figure  12. 

For  some  of  the  test  runs,  the  digital  timer  located  immediately 
below  the  control  panel,  as  shown  in  Figure  12,  was  used  in  conjunction 
with  the  oscillograph  and  the  oscilloscope  to  determine  the  reaction 
propagation  velocity.  This  device  was  sc  designed  as  to  utilize 
the  thermocouple  analog  output  voltage  signals  to  start  and  stop  a 
series  of  digital  counters  which  measured  and  displayed  the  time 
elapsed  between  the  arrival  of  the  reaction  front  at  subsequent 
thermocouple  locations.  Electrical  noise  problems,  however,  rendered 


Figure  12  Pho 
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the  operation  of  the  digital  timer  highly  unreliable,  and  the  use  of 
this  device  was  limited  to  the  occasional  corroboration  of  the  records 
produced  by  the  oscillograph  and  the  oscilloscope  The  vacuum 
pump  located  below  the  digital  timer,  r.nd  the  argon  tank  positioned 
behind  the  control  panel  are  also  visible  in  Figure  12.  The  igniter 
was  energized  by  a 12  volt  output,  series-parallel  arrangement  of 
four,  6-volt  "Hot  Shot"  batteries  also  located  behind  the  control 
panel. 

In  some  of  the  test  runs,  it  was  attempted  to  measure  and 
record  the  pressure  occurring  t 'n  the  sample  tube  while  the 
reaction  was  taking  place.  .er  to  accomplish  this,  two  of  the 

304  stainless  steel  tie  rode  in  the  test  apparatus  were  replaced  by 
aluminum  tie  rods  incorporating  reduced  cross-sectional  area  sections 
which  were  strain  gauged.  The  two  remaining  stainless  steel  tie 
rods  were  allowed  to  remain  in  place  but  were  installed  loosely, 
acting  as  overpressure  protection  only.  The  strain  gauges  on  the  two 
aluminum  tie  rods  were  wired  in  series  in  order  to  minimize  any 
unsymmetrical  effects,  and  were  included  as  one  arm  of  a standard 
Wheatstone  bridge.  The  bridge  circuit  was  powered  by  a 6 volt  dry 
cell  battery  and  its  output  was  suitably  amplified  and  recorded  on 
che  CEC  oscillograph.  The  system  was  calibrated  prior  to  each  test 
run  by  applying  known  pressures  to  the  interior  of  ar.  empty  sample 
tube  and  observing  the  resulting  output  traces  on  the  oscillograph 


record . 
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3 . 2 Experimental  Procedure 

The  preparation  for,  and  the  conduction  of  a typical  test 
run  took  place  in  the  following  manner:  An  empty  sample  tube,  three 

thermocouple  probes,  and  a core  rod  are  shown  in  the  photograph  in 
Figure  13.  The  core  rod,  having  been  given  a light  coating  of 
mineral  oil  to  prevent  lithium  from  adhering  to  its  surface,  was 
first  installed  in  the  empty  sample  tube  by  screwing  the  nut  at  che 
end  of  the  rod  into  the  igniter  hole  in  the  igniter  plate.  The  core 
rod  was  thus  positioned  along  the  center  of  the  sample  tube.  The 
thermocouple  probes  were  then  inserted  through  the  holes  in  the  wall 
of  the  sample  tube  until  their  measuring  junctions  contacted  the 
surface  of  the  core  rod,  and  were  secured  in  this  position  with  high 
temperature  epoxy.  This  sample  tube-thermccouple  probe-core  rod 
assembly  was  then  placed  in  an  inert-atmosphere  glove  box  (Vacuum 
Atmospheres  Model  HE-133-5)  along  with  the  cap  and  shutoff  valve,  the 
retainer,  and  the  tie  rods,  Under  an  atmosphere  of  dry  argon, 
liquid  lithium  at  approximately  480  K was  then  poured  into  the  empty 
sample  tube,  surrounding  the  core  rod  and  the  thermocouple  probes. 
After  the  lithium  had  solidified  and  '"eturned  to  roo-  temperature , 
the  core  rod  was  t emoved  and  a dummy  igniter  installed  in  the  igniter 
plate.  The  PTFE  sample  was  then  installed,  followed  by  the  small 
quantity  of  lithium  fluoride  The  wafer  having  been  positioned  in  the 
end  of  the  sample  tube,  the  cap,  the  retainer,  and  the  tie  rode  were 
put  in  place.  The  shutoff  valve  war.  then  closed,  and  the  test 
apparatus  was  removed  from  the  glove  box 

The  test  apparatus  was  then  secured  to  the  lid  of  the 


safety  vessel,  the  thermocouple  leads  were  connected  to  the  thermocouple 
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jack  panel  through  the  extension  wires,  and  the  solenoid  valve  and 
attendant  plumbing  were  installed.  The  dummy  igniter  was  replaced 
by  a live  igniter,  the  igniter  power  lead  was  connected,  and  the 
shuto£f  valve  was  opened.  The  teat  apparatus  was  then  lowered  into 
the  safety  vessel  by  means  of  an  electric  hoist,  and  the  lid  of  the 
vessel  was  secured. 

The  solenoid  valve  was  then  opened  and  the  sample  tube  and 
the  irterlor  of  che  safety  vessel  were  evacuated.  After  the 
attainment  of  the  desired  pressure  level  within  the  sample  tube  and 
the  safety  vessel,  the  solenoid  valve  was  closed,  as  were  the  shutoff 
valves  on  the  lid  of  the  vessel  and  those  adjacent  to  the  pressure 
gauges  on  the  control  panel  The  oscillograph  chart  drive  was  then 
started  and  the  igniter  was  fired,  electric  power  to  the  igniter 
being  also  utilized  to  trigger  the  oscilloscope  sweep 

A sketch  of  the  temperature-time  record  produced  by  the 
oscillograph  tor  a typical  test  run  is  shown  in  Figure  14-  As  seen 
in  this  sketch,  che  arrival  of  the  reaction  front  at  each  thermocouple 
location  was  attended  by  a rapid  temperature  rise,  the  reaction  front 
reaching  the  location  of  the  thermocouple  nearest  the  igniter 
approximately  10  ms  after  electric  powe;  was  applied  to  the  igniter. 
From  timing  traces  present  cn  the  oscillograph  record,  the  time 
Intervals  between  the  arrival  of  the  reaction  front  at  each  of  the 
thermocouple  locations  were  ascertained,  and  the  reaction  propagation 
velocity  was  determined  by  dividing  the  distance  between  the 
thermocouple  measuring  junctions  by  the  time  required  for  the  reaction 
front  to  travel  from  one  thermocouple  location  to  the  next-  The 


total  time  required  for  the  reaction  front  to  travel  from  the 


TIME *- 

Figure  14  Sketch  of  Temperature-Time  Record  for  a Typical  Test 
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thermocouple  nearest  the  igniter  to  the  thermocouple  nearest  the 
cap  was  determined  in  a similar  manner  from  the  screen  of  the  storage 
oscilloscope . 

When  the  excess  lithium  and  the  reaction  products  contained 
in  the  sample  tube  had  cooled  to  a reasonable  temperature,  the  sample 
tube  and  the  safety  vessel  were  backfilled  to  atmospheric  pressure 
with  argon  gas.  The  lid  of  the  safety  vessel  and  the  attached  test 
apparatus  were  then  removed  from  the  vessel.  Finally,  the  test 
apparatus  was  disassembled,  the  sample  tube  and  its  contents  being 
discarded,  ^nd  the  cap,  the  shutoff  valve,  and  the  associated 
plumbing  be.lra  cleaned  of  residual  lithium  by  water  immersion  in 
preparation  fcr  a future  test  run. 

It  should  be  pointed  out,  with  reference  to  Figure  14, 
that  the  accuracy  with  which  the  reaction  propagation  velocity  could 
be  measured  was  seriously  compromised  by  the  nature  of  the  temperature- 
time record.  In  order  that  the  time  intervals  between  the  arrival 
of  the  reaction  front  at  subsequent  thermocouple  locations  could  be 
determined,  it  was  necessary  t:  measure  very  small  physical  distances 
between  the  rapid  temperature  rises  displayed  on  the  record,  which 
introduced  an  estimated  +15%  uncertainty  in  the  determination  of  the 
reaction  propagation  velocity  Since  it  was  not  possible  to  further 
increase  the  oscillograph  chart  speed,  which  would  have  allowed 
better  resolution  in  these  measurements,  the  developr.c.i;.  cl  Jie 
previously  described  digital  timer  was  undertaken,  it  being 
theoretically  possible  to  greatly  increase  the  accuracy  oi  the 
measurements  through  the  use  of  this  device  Although  much  effort 
was  expended  in  attempts  to  rectify  the  situation,  the  digital  timer 
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consistently  malfunctioned  under  actual  test  conditions.  For  this 
reason,  the  temperature-time  records  produced  by  the  oscillograph 
were  used  to  determine  the  reaction  propagation  velocities  in  spite 


of  the  inherent  inaccuracies  present. 


CHAPTER  IV 

EXPERIMENTAL  AND  THEORETICAL  RESULTS 

4.1  Experimental  Conditions 

In  the  experimental  phase  of  this  study,  the  effects  of 
varying  the  diameter  and  density  of  the  PTFE  sample  on  the  lithium- 
PTFE  reaction  propagation  velocity  were  investigated.  Density 
variation  was  accomplished  by  using  solid  PTFE  rod  for  the  highest 
density  condition,  and  PTFE  resin  at  several  degrees  of  compaction 
for  the  less  dense  situations.  The  use  of  core  rods  of  different 
sizes  in  the  test  preparation  process  allowed  the  diameter  of  the 
PTFE  sample  to  be  changed.  Test  runs  were  conducted  with  solid 
PTFE  rod  samples  of  7.94  nun  (5/16  inch),  6.35  mm  (1/4  inch),  4.76 
mm  (3/16  inch),  and  3.18  mm  (1/8  inch)  diameters,  and  with  PTFE  resin 
samples  having  diameters  of  9.52  mm  (3/8  inch),  and  6.35  mm.  Lithium 
metal  obtained  from  the  Lithium  Corporation  of  America,  with  99.8% 
minimum  guaranteed  purity,  was  used  in  all  test  runs. 

The  PTFE  resin  utilized  in  this  investigation  was  Teflon 
Fluorocarbon  resin  manufactured  by  E J duPont  de  Nemours  and 
Company.  Teflon  8A  resin,  having  an  average  particle  size  of  0.52 
mm,  vas  used  in  the  majority  of  the  PTFE  resin  test  runs.  Teflon  8 
resin,  and  Teflon  9B  resin,  with  average  particle  sizes  of  0.65  mm 
and  0.55  mm,  respectively,  were  used  in  a small  number  of  test  runs 
in  order  to  assess  particle  size  effects.  The  largest  portion  of  the 
Teflon  8A  resin  was  obtained  from  duPont  through  normal  commercial 
channels,  while  the  remaining  Teflon  8A  resin,  as  well  as  the  Teflon  8 
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and  9B  resins,  was  furnished  by  the  Technical  Service  Laboratory, 
Plastics  Department,  Fluorocarbons  Division,  at  duPont.  The 
apparent  bulk  density  of  these  PTFE  resins  in  their  uncompacted 
form  was  of  the  order  of  650  kg/m  . As  a result  ot  the  sample 
preparation  process  for  the  PTFE  resin  test  runs,  which  involved 
lightly  tamping  the  resin  into  the  hole  in  the  lithium  sample  left 
by  the  removal  of  the  core  rod,  the  density  of  the  PTFc.  resin  samples 
used  in  these  test  runs  was  increased  to  values  m the  range  of 
730-1200  kg/m3.  In  order  to  preclude  the  possibility  of  a premature 
ignition  caused  by  excessive  lithium-PTFE  contact  pressure,  no 
attempts  were  made  to  pack  the  PTFE  resin  tc  densities  greater  than 
1200  kg/m3. 

In  approximately  SOX  of  the  test  runs  involving  solid 
PTFE  rod,  general  purpose  grade  extruded  PTFE  rod  obtained  commerically 
from  the  Cadillac  Plastic  and  Chemical  Company  was  utilized.  In  the 
remainder  of  these  test  runs,  extruded  PTFE  rod  furnished  by  duPont' s 
Technical  Service  Laboratory  was  used-  The  rod  obtained  from  duPont 
consisted  ot  Semples  fabricated  from  Teflon  8,  8A  and  9B  resins, 
these  resins  being  from  the  same  lots  as  were  the  resins  furnished 
by  auPont  used  in  the  PTFE  resin  test  runs.  The  PTFE  rod  samples, 
as  used  in  all  but  one  of  the  test  runs,  were  machine-turned  to  a 
diameter  nominally  0.25  mm  (0.010  inch)  smaller  than  the  core  rod 
diameter  in  order  to  avoid  interference  with  the  lithium  during 
installation . 

In  all  but  one  case,  the  lithium-PTFE  reaction  test  runs 
conducted  in  this  investigation  were  carried  out  with  the  reaction 
sample  under  an  initial  vacuum.  The  ambient  pressure  level  within 
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the  3ample  tube  immediately  prior  to  the  firing  of  the  igniter  in 

-3 

these  test  runs  was  approximately  10  bar. 

4 2 Test  Runs  with  9.52  nan  Diameter  PTFE  Samples 

The  test  runs  conducted  with  9.51  mm  f 3 8 inch)  diameter 

PTFE  samples  in  all  cases  involved  the  use  of  PTFE  resin  These  test 

runs  are  summarized  in  Table  1.  The  first  5 runs  listed  in  Table  1, 

Runs  W-l  through  W-5,  are  the  preliminary  tests  described  in  Chapter 

1 For  only  these  five  runs,  the  reaction  sample  preparation  process 

was  not  carried  out  in  the  inert-atmosphere  glove  box,  but  rather  was 

done  in  the  ambient  atmosphere,  the  sample  being  continuously  bathed 

in  argon  gas  in  order  to  minimize  air  contamination  of  the  lithium. 

Further,  no  lithium  fluoride  was  present  in  the  reaction  sample,  and 

the  wafer  at  the  end  of  the  sample  tube  wss  not  utilized. 

Commercially  obtained  Teilon  8A  resin  at  densities  in  the  range  780- 
3 

900  kg/m  was  used  in  these  test  runs,  yielding  PTFE-lithium  mass 
ratios  from  0 226  to  0261.  The  reaction  propagation  velocities 
listed  in  Table  1 for  these  five  test  runs,  as  well  as  for  a<.l  other 
runs  with  9. 52  mm  diameter  PTFE  samples,  were  determined  by  measuring 
the  time  required  for  the  reaction  front  to  travel  from  the 
thermocouple  located  nearest  the  igniter  tc  the  thermocouple  nearest 
the  cap.  For  sll  runs  with  S-52  nun  diameter  PTFE  samples,  the  reaction 
propagation  velocities  measured  over  the  pcrtion  ci  '.he  sample 
between  the  thermocouple  nearest  the  igniter  and  the  thermocouple  at 
the  midpoint  of  the  sample,  and  between  this  thermocouple  and  the 
t Her mo couple  nearest  the  tap  were  not  significantly  different  from 
the  reaction  propagation  velocity  determined  by  using  only  the 
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thermocouple  nearest  the  igniter  and  the  thermocouple  neatest  the  cap. 
Under  these  relatively  constant  velocity  conditions,  the  reaction 
propagation  velocity  measurements  made  using  the  two  thermocouples 
nearest  the  ends  of  the  reaction  sample  were  considered  more  reliable 
due  to  the  smaller  effect  of  errors  in  the  axial  placement  of  the 
thermocouples.  In  these  five  test  runs,  the  reaction  pr 'pagation 
velocities  were  determined  through  the  use  of  the  CEC  Oscillograph 
alone.  As  seen  from  Table  1,  four  of  these  five  preliminary  test 
runs  were  successful.  In  these  runs,  the  reaction  front  propagated 
over  the  entire  length  of  the  reaction  sample,  alJ owing  the  reaction 
propagation  velocity  to  be  determined.  During  Run  W-3 , however,  the 
reaction  front  propagated  through  only  a portion  cf  the  sample  before 
propagation  ceased.  This  phenomenon  was  attributed  to  nonuniform 
packing  of  the  PTFE  resin,  as  previously  described 

In  the  next  five  test  runs  with  the  9 . 52  mm  diameter  PTFE 
resin  sample.  Runs  M-l  through  M-5  in  Table  1,  it  was  attempted  to 
measure  the  reaction  propagation  velocity  by  means  of  the  digital 
timer  alcne , For  these  runs,  commercially  obtained  Teflon  8A  resin 
was  used  at  densities  from  830  to  1000  kg/m^  The  PTFE-lichium 
mass  ratios  present  in  these  runs,  as  well  as  in  all  ether  runs  in 
which  lithium  fluoride  was  used  to  inhibit  the  flow  :i  high 
temperature  lithium  into  the  cap  and  associatei  plumbing,  were 
computed  by  neglecting  that  portion  of  the  lithium  sample  which 
surrojnded  the  lithium  fluoride  The  PTFE-lithium  mass  ratios  fo' 
runs  M-l  through  M-5  were  thus  0 242  to  0.295.  In  each  of  these  test 
runs,  the  digital  time  maltunc tioued , and  no  reaction  propagation 
velocity  data  was  obtained.  In  ail  subsequent  test  runs  listed  in 
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Table  1,  as  well  as  all  test  runs  with  PTFE  sample  diameters  other 
than  9.52  mm,  the  CEC  Oscillograph  was  used  as  the  primary  means  of 
measuring  the  reaction  propagation  velocity,  with  the  storage 
oscilloscope  being  utilized  for  corroboration  of  the  results. 

In  an  effort  to  make  the  results  of  the  test  runs  less 
sensitive  to  suspected  minor  variations  in  the  surface  finish  of  the 
lithium  at  the  lithium-PTFE  interface,  for  Runs  C-l  through  C-5 
the  thermocouple  measuring  junctions  were  moved  1.6  mm  (1/16  inch) 
away  from  the  lithium-PTFE  interface,  into  the  lithium  sample. 
Commercially  obtained  Teflon  8A  resin  at  densities  between  990  and 

3 

1080  kg/m  was  used  in  these  five  runs,  the  resulting  PTFE-lithium 
mass  ratios  being  0.278  to  0.301.  In  three  of  these  test  runs, 

C-2,  C-3  and  C-4  as  shown  in  Table  1,  the  reaction  front  propagated 
over  the  entire  length  of  the  reaction  sample.  During  Run  C-l  the 
propagation  of  the  reaction  ceased  before  the  end  of  the  sample  had 
been  reached,  and  during  Run  C-5,  the  reaction  was  not  initiated  by 
the  firing  of  the  igniter.  Of  the  three  test  runs  in  which  the 
propagation  was  complete,  only  Run  C-2  yielded  reaction  propagation 
velocity  data.  The  analog  records  from  the  oscillograph  for  the 
remaining  two  test  runs  indicated  that  the  reaction  front  arrived  at 
the  thermocouple  located  nearest  the  cap  before  arriving  at  the 
thermocouple  located  nearest  the  igniter.  The  atmosphere  present 
within  the  sample  tube  immediately  prior  to  the  firing  of  the  igniter 
m Run  C-4  was  argon  at  1 bar  pressure,  as  compared  to  the  10  bar 
vacuum  present  in  all  other  runs.  The  nonsequential  heating  of 
the  thermocouples  was,  however,  not  effected  by  this  change  in 
initial  test  conditions.  It  being  physically  impossible  for  the 
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reaction  front  to  arrive  at  the  thermocouple  farthest  from  the 
igniter  before  arriving  at  the  thermocouple  nearest  the  igniter, 
assuming  the  thermocouples  were  located  at  the  same  radial  position, 
the  results  of  these  test  runs  indicated  that  the  reaction  front  is 
very  thin  in  the  radial  direction,  and  that  small  errors  in  the 
positioning  of  the  thermocouple  measuring  junctions  at  some  distance 
from  the  lithium-PTFE  interface  are  very  critical.  For  this  reason, 
the  reaction  propagation  velocity  measurement  from  Run  C-2  was  Judged 
not  reliable,  and  the  thermocouple  measuring  junctions  were  located 
at  the  llthium-PTFE  inteiface  for  all  remaining  test  runs  summarized 
in  Table  1,  as  well  as  for  all  test  runs  with  PTFE  sample  diameters 
other  than  9.52  mw. 

Test  Runs  C-6  and  C-7  were  conducted  with  commercially 

3 

obtained  T»  :lor.  8A  resin  at  1030  kg/m  density,  and  PTFE-lithium  mass 
ratios  of  v.286  an.'.  0.285,  respectively,  As  seer,  from  Table  1,  the 
quantities  of  lxthium,  PTFE,  and  lithium  fluoride  used  in  these  two 
runs  were  virtually  Identical.  The  reaction  propagation  velocity 
measured  in  Run  C-6  vas  3?  m/s,  and  that  in  Run  C-7,  36  m/s  Fr.m 
tnese  result*;,  ic  is  indicated  that  under  similar  experimental 
conditions,  the  reaction  propagation  measurements  are  quite  repeatable. 

In  the  next  cwo  test  runs  shown  in  Table  3 , an  attempt  was 
made  to  assess  the  effect  cf  air  contaminaticn  c f the  lithium  sample 
on  the  reaction  propagation  velocity.  These  tests.  Runs  C-8  and 
C-9,  resembled  Runs  C-6  and  C-7  in  that  commercially  obtained  Teflon 
8A  resin  was  again  used,  and  the  amounts  of  lithium,  PTFE,  and  lithium 
fluoride,  the  PTFE  resin  densities,  and  the  PTFE-lithium  mass  ratios 
were  similar-  After  the  core  rod  had  been  withdrawn,  but  before  the 
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installation  of  the  PTFE  resin  during  the  preparation  for  Run  C-8, 
the  lithium  sample  in  the  sample  tube  was,  however,  removed  from  the 
inert-atmosphere  glove  box  and  exposed  to  the  ambient  air  for  a period 
of  10  seconds.  The  sample  tube  was  then  replaced  in  the  glove  box 
and  the  sample  preparation  process  was  continued  In  the  normal  manner. 
The  reaction  propagation  velocity  measured  in  this  test  run  was 
32  m/s.  The  preparation  for  Run  C-9  was  carried  out  in  a like  fashion, 
the  lithium  sample  in  this  case  being  exposed  to  the  ambient  air  for 
a period  of  60  seconds,  and  the  reaction  propagation  velocity  being 
55  m/s.  The  reaction  propagation  velocity  in  Run  C-8,  32  m/s, 
compares  favorably  with  that  observed  in  the  previous  two  similar 
test  runs,  while  the  reaction  propagation  velccity  in  Run  C-9  was 
considerably  higher.  It  is  therefore  indicated  that  the  contamination 
of  the  lithium  surface  caused  by  prolonged  air  exposure  has  a 
noticeable  effect  on  the  reaction  propagation  velocity  The  exact 
mechanirm  through  which  the  reaction  propagation  process  is  altered 
in  this  situation  is  open  to  speculation  It,  howe”er,  seems 
likely  that  the  presence  of  the  layer  of  lithium  hydroxide  iormed  on 
the  surface  of  the  lithium  as  a result  of  its  expcture  to  moist 
ambient  air  (3)  might  well  cause  a significant  change  in  the  thermal 
conductivity  of  the  lithium  at  the  lrthi um-PTFE  interlace,  and  that 
this  difference  ir.  thermal  conductivity  might  be  reflected  by  a change 
in  rho  reaction  propagation  velocity. 

Test  Run  C-10  was  conducced  in  :rder  that  the  effect  of 
the  use  of  Teflon  8A  resin  obtained  from  a different  scarce  might  be 
assessed.  The  resin  utilized  in  this  run  was  furnished  by  the 
Technical  Service  Laboratory  at  duPont,  and  was  capable  of  being 
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easily  camped  to  higher  density  than  was  the  Te£lon  8A  resin 
obtained  from  duPont  commercially.  As  shown  in  Table  1,  the  resin 
density  was  1130  kg/m\  tha  PTFE-lithium  mass  ratio  was  0.312,  and 
the  reaction  propagation  velocity  was  51  m/s  for  this  test  run.  Since 
the  resin  density  and  the  PTFE-lithium  mass  ratio  for  Run  C-10  were 
approximately  10%  higher  than  in  Runs  C-6  and  C-7,  made  with 
commercially  obtained  Teflon  8A  resin,  it  is  impossible  to  ascertain 
whether  the  increase  in  the  reaction  propagation  velocity  was  due  to 
the  change  in  test  conditions  alone,  or  to  some  difference  in  the 
PTFE  resin  itself. 

In  Run  C-10,  the  pressure  level  present  within  the  sample 
tube  while  the  reaction  was  caking  place  was  measured  by  means  of 
the  strain-gauged  aluminum  tie  rods  described  in  Chapter  III,  The 
maximum  pressure  level  was  12.4  bar  (180  psia) . This  pressure 
maximum  occurred  at  approximately  the  same  time  that  the  reaction 
front  reached  the  measuring  junction  of  the  thermocouple  located 
nearest  the  igniter,  indicating  that  the  pressure  wave  propagated 
through  the  reaction  sample  more  rapidly  than  did  the  reaction 
front . The  pressure  '.'ithin  the  sample  tube  decreased  approximately 
10%  from  tire  maximum  during  the  time  required  for  the  reaction  front 
to  travel  from  tire  thermocouple  located  nearest  the  igniter  to  the 
thermocouple  located  nearest  the  cap. 

The  final  test  run  with  the  9.5?.  mm  diameter  PTFE  sample, 

Run  C-J1 , was  carried  out  in  order  to  evaluate  cna  effect  of  the  uc-e 
rt  a different  type  PTFE  resin  on  the  reaction  propagation  velocity. 

As  shown  in  Table  Teflon  * resin  obtained  from  duPont 's  Technical 
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Service  Laboratory  was  utilized  in  this  run,  the  resin  density  and  the 
PTFE-lithium  mass  ratio  being  again  higher  than  in  previous  test  runs. 
In  Run  C-ll,  the  reaction  was  not  initiated  by  the  firing  of  the 
igniter.  Whether  this  phenomenon  was  caused  by  a faulty  igniter, 
or  some  property  of  the  Teflon  8 resin  is  unknown. 

Of  the  total  of  21  test  runs  conducted  with  the  the  9.52  mm 
diameter  PTFE  resin  sample  and  summarized  in  Table  1,  then,  reaction 
propagation  velocity  data  was  obtained  in  10  runs,  and  9 of  these 
10  produced  data  judged  reliable.  The  successful  test  runs  3n  which 
commercially  obtained  Teflon  8A  resin  was  utilized  included  4 
preliminary  test  runs,  2 air  exposure  runs,  and  2 "baseline"  runs 
with  nearly  identical  experimental  conditions.  In  these  baseline 
runs.  Runs  C-6  and  C-7,  the  average  measured  reaction  propagation 
velocity  was  36.5  m/s  at  an  average  PTFE-lithium  mass  ratio  of 
0.2853.  The  results  of  the  short  duration  air  exposure  run  (Run  C-8) 
are  in  reasonable  agreement  with  those  of  the  baseline  runs,  while 
the  results  of  the  longer  duration  air  exposure  run  (Run  C-9)  are 
not.  The  preliminary  test  runs  were  conducted  at  an  average  PTFE- 
lithium  mass  ratio  of  0.245,  and  yielded  an  average  measured  reaction 
propagation  velocity  of  36  m/s.  Although  the  actual  exposure  of  the 
lithium  to  the  ambient  air  during  the  preparation  for  these 
preliminary  runs  was  probably  of  short  duration,  the  results  of  these 
runs  and  their  agreement  with  those  of  the  baseline  runs,  are 
nevertheless  of  questionable  reliability. 

The  reaction  propagation  velocity  in  the  test  run  in  which 
Teflon  8A  resin  obtained  from  the  Technical  Service  Laboratory  at 
dcPont  was  used  fRun  C-10)  was  approximately  40%  higher  then  that 
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observed  in  the  baseline  runs;  the  exact  cause  of  this  increase  is 
not  known.  In  this  run,  pressure  measurements  indicated  that,  at 
least  in  this  particular  situation,  the  velocity  with  which  the 
pressure  wave  moves  through  the  sample  exceeds  chat  with  which  the 
reaction  front  moves.  The  reaction  was  unable  to  be  initiated  by  normal 
means  in  the  one  run  conducted  with  Teflon  8 resin,  although  no 
obvious  explanation  for  this  phenomenon  exists. 

4 . 3 Test  Runs  with  6.35  mm  Diameter  FIFE  Samples 

In  the  test  runs  conducted  with  6.35  mm  (1/4  inch)  diameter 
PTFE  samples,  PTFE  in  the  form  of  resin,  as  well  as  in  the  form  of 
solid  rod  was  utilized.  These  test  runs  are  summarized  in  Table  2 and 
Table  3,  respectively.  The  reaction  sample  preparation  process  for 
all  runs  with  6.35  nun  diameter  PTFE  samples  was  carried  out  in  the 
inert -atmosphere  glove  box,  and  the  analog  output  records  from  the 
CEC  Oscillograph  were  used  as  the  primary  means  of  determining  the 
reaction  propagation  velocity,  data  from  the  storage  oscilloscope 
being  used  for  corroboration.  As  was  the  case  in  the  test  runs  with 
the  9,52  mm  diameter  PTFE  samples,  the  reaction  propagation  velocity 
in  all  runs  with  6.35  mm  diameter  PTFE  samples  was  computed  by  using 
data  obtained  from  the  two  thermocouples  nearest  the  ends  of  the 
reaction  sample,  it  being  again  observed  that  the  reaction  propagation 
velocity  was  reasonably  constant  over  the  entire  instrumented  length 
of  these  samples.  The  thermocouple  measuring  junctions  for  these 
runs  were  located  at  the  lithium-PTFE  interface,  the  PTFE-lithium 
mass  ratios  were  computed  by  neglecting  the  lithium  surrounding  the 
lithium  fluoride  adjacent  to  the  cap,  and  the  pressure  levels  within 
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the  sample  tubes  Immediately  prior  to  the  firing  of  the  igniter  were 

_3 

approximately  10  bar. 

As  shown  in  Table  2,  all  test  runs  conducted  with  PTFE 
resin  samples  of  6.35  mm  diameter  involved  the  use  of  Teflon  8A 
resin.  In  the  first  two  runs  listed  in  Table  2,  Runs  D-l  and  D-2, 
commercially  obtained  Teflon  8A  resin  was  utilized  at  densities  of 
1060  and  1050  kg/m  , respectively.  The  amounts  of  lithium,  PTFE, 
and  lithium  fluoride  present  in  these  two  runs  were  similar, 
as  were  the  PTFE-lithium  mass  ratios.  The  reaction  propagation 
velocity  measured  in  Run  D-l  was  42  m/s,  and  that  in  Run  D-2,  62  m/s. 

The  next  two  test  runs  with  the  6.35  mm  diameter  PTFE  resin 

sample  were  carried  out  in  order  to  investigate  the  effect  of  reduced 

PTFE  resin  density  on  the  reaction  propagation  velocity.  As  shown 

in  Table  2,  test  runs  D-3  and  D-4  involved  equal  amounts  of  lithium, 

and  similar  amounts  of  lithium  fluoride.  These  runs  were  conducted 

3 

with  PTFE  resin  densities  of  730  and  830  kg/m  , and  PTFE-lithium  mass 
ratios  of  0.083  and  0.096 , respectively , Teflon  8A  resin  obtained  from 
the  Technical  Service  Laboratory  at  duPont  being  utilized.  Since  the 
thermocouple  measuring  junctions  were  heated  nonsequentially  in 
Runs  D-3  and  D-4,  the  thermocouple  nearest  the  cap  indicating  a 
sharp  temperature  rise  before  the  thermocouple  nearest  the  igniter, 
no  reaction  propagation  velocity  data  resulted.  Post-test  run 
inspection  of  the  reaction  sample  showed  that  the  PTFE  resin  had  been 
blown  toward  the  opposite  end  of  the  sample  tube  by  the  firing  of  the 
igniter.  The  exact  resin  density  below  which  the  reaction  will  not 
propagate  normally  is,  no  doubt,  a function  of  the  type  of  PTFE  resin 


1 

} 

i 


\ 

i 


utilized.  However,  in  these  two  test  runs  it  was  demonstrated  that, 
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at  sufficiently  low  PTFE  resin  densities,  the  configuration  of  the 
PTFE  sample  is  altered  by  the  action  of  the  igniter  to  such  an  extent 
that  normal  reaction  propagation  is  prohibited. 

Test  Runs  D-5  and  D-6  were  the  final  two  runs  involving 

6.35  mm  diameter  PTFE  samples,  as  shown  in  Table  2.  Commercially 

obtained  Teflon  8A  resin  was  used  in  these  runs,  and  experimental 

conditions  were  made  similar  to  those  present  in  Runs  D-l  and  D-2 

in  an  attempt  to  establish  a more  meaningful  baseline  for  the  PTFE 

resin  test  runs  in  this  PTFE  sample  diameter.  In  Run  D-5,  the  reaction 

was  not  initiated  by  the  firing  of  the  igniter.  This  phenomenon 

was  attributed  to  the  inadvertent  exposure  of  the  lithium  sample  to 

ambient  air  for  a period  of  several  hours  during  the  sample 

preparation  process.  Test  Run  D-6,  in  wnich  the  PTFE  resin  density 
3 

was  1080  kg/m  , and  the  PTFE  mass  ratio  was  0.123,  yielded  a measured 
reaction  propagation  velocity  of  49  m/s.  The  pressure  level  within 
the  sample  tube  during  the  reaction  propagation  process  was  also 
measured  in  this  run,  the  maximum  pressure  present  being  15  9 bar 
(230  psia).  In  Run  D-6,  no  pressure  was  evidenced  until  after  the 
reaction  front  had  passed  the  measuring  junction  of  the  thermocouple 
located  nearest  the  cap,  indicating  that,  in  this  run,  the  reaction 
front  propagated  through  the  sample  at  least  as  rapidly  as  did  the 
pressure  wave. 

From  the  results  of  the  test  runs  conducted  with  6 . 35  mm 
diameter  PTFE  resin  samples,  as  shown  in  Table  2,  it  was  observed 
that  the  reaction  front  will  not  propagate  normally  at  sufficiently 
low  PTFE  resin  densities.  The  three  successful  baseline  runs  (Runs 
D-l,  D-2  and  D-6),  with  an  average  PTFE-iithium  mass  ratio  ci  0.121, 
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yielded  an  average  measured  reaction  propagation  velocity  of  51  m/s, 
and  pressure  measurements  made  during  Run  D-6  indicated  that  under 
the  experimental  conditions  shown  in  Table  2,  the  reaction  front  moves 
through  che  reaction  sample  at  a velocity  net  less  than  that  of  the 
pressure  wave. 

Shown  in  Table  3 is  a summary  of  the  test  runs  conducted 
with  6.35  mm  diameter  solid  PTFE  rod  samples.  For  these  runs, 

i 

commercial  PTFE  rod  obtained  from  the  Cadillac  Plastic  and  Chemical 

Company,  as  well  as  PTFE  rods  fabricated  from  Teflon  8,  Teflon  8A, 

and  Teflon  9B  resins  by  the  Technical  Service  Laboratory  at  duPont, 

were  used  in  the  construction  of  the  PTFE  samples.  As  in  the  test 

runs  with  PTFE  resin  samples  of  the  same  diameter,  the  pressure  level 

within  the  sample  tube  at  the  time  of  the  firing  of  the  igniter  in 

_3 

these  runs  was  approximately  10  bar.  The  reaction  samples  were 
again  prepared  in  the  inert-atmosphere  glove  box,  che  thermocouple 
measuring  junctions  being  placed  at  the  lithium-PTFE  interface.  The 
reaction  propagation  velocities  in  these  runs  were  determined  from 
the  analog  output  records  from  the  CEC  Oscillograph  and  were  confirmed 
by  the  data  from  the  storage  oscilloscope,  and  the  lithium  surrounding 
the  LiF  at  the  end  of  the  sample  was  neglected  in  the  computation 
of  the  PTFE-lithium  mass  ratios. 

Test  Runs  E-i,  E-2,  and  E-8  were  conducted  with  6.35  mm 
diameter  commercial  solid  PTFE  rod  samples,  as  shown  in  Table  3.  For 
Run  E-l,  the  PTFE  sample  was  machined  tc  e diameter  0,50  mm  (0-020 
inch)  smaller  than  the  core  rod  diamete'  in  order  to  insure  that  no 
interference  between  the  lithium  and  the  PTFE  samples  during  the 
installation  of  the  PTFE  would  occur-  A rather  large  gap  between  the 
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lithium  and  the  PTFE  at  the  lithium-PTFE  interface  was  thus  created. 
Although  the  reaction  front  propagated  in  a normal  manner  in  this 
test  run,  at  a measured  reaction  propagation  velocity  of  70  m/s,  this 
result  was  judged  unreliable  due  to  the  excessive  clearance  between 
the  lithium  and  PTFE  samples.  The  PTFE  samples  for  all  other  test 
runs  involving  solid  PTFE  rod  were  machined  to  a diameter  nominally 
0.25  mm  (0.010  inch)  smaller  than  the  core  rod,  which  reduced  the 
clearance  between  the  samples  while  still  preventing  an  interference 
fit.  Test  Runs  E-2  and  E-8  were  conducted  with  PTFE-lithium  mass 
ratios  of  0.224  and  0.219,  respectively.  The  measured  reaction 
propagation  velocity  in  Run  E-2  was  57  m/s,  and  that  in  Run  E-8, 

41  m/s.  From  a leak  which  developed  during  Run  E-8,  a quantity  of 
hot  lithium  was  deposited  in  the  vicinity  of  one  of  the  strain-gauged 
aluminum  tie  reds.  The  thermocouple  extension  wires  were  uneffected, 
however,  and  the  reaction  propagation  velocity  measurement  was  judged 
reliable . 

In  Runs  E-3,  E-4  and  E-5,  the  PTFE  samples  were  constructed 
from  solid  rods  fabricated  from  Teflon  8A  resin  by  the  Technical 
Service  Laboratory  at  duPont-  During  Run  E-4,  a lithium  leak  onto 
the  thermocouple  extension  wires  prevented  reaction  propagation  data 
from  being  obtained.  As  a result  of  this  lithium  leak,  and  the 
attendant  destruction  of  the  thermocouple  extension  wires,  it  was 
not  possible  to  ascertain  whether  the  propagation  of  the  reaction 
front  occurred  normally.  In  addition,  it  was  later  determined  that 
cutting  oil  had  been  used  on  the  PTFE  sample  during  final  sizing, 
it  being  felt  that  the  presence  of  a small  quantity  or  residual 
ell  on  the  PTFE  sample  might  well  effect  tte  reaction  propagation 
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process.  In  Runs  E-3  and  E-5,  the  reaction  front  propagated  through 
the  sample  in  a normal  manner  at  measured  propagation  velocities  of 
42  m/s  and  50  m/s .respectively , as  shown  in  Table  3.  For  Run  E-5, 
two  additional  thermocouples  were  placed  at  the  same  axial  position 
as  the  thermocouple  at  the  midpoint  of  the  sample  tube,  but  at 
different  circumferential  locations  around  the  lithium-PTFE  interface. 

The  reaction  front  was  observed  to  arrive  at  the  measuring  junctions 

\ 

of  these  three  thermocouples  simultaneously,  indicating  that  the 
propagation  of  the  reaction  front  is  axisymmetric, 

Solid  PTFE  rod  fabricated  from  Teflon  8 resin  at  duPont's 
Technical  Service  Laboratory  was  used  in  the  construction  of  the 
PTFE  samples  for  Runs  E-6  and  E-9 , as  shown  in  Table  3,  In  both  of 
these  runs,  the  reaction  front  did  not  propagate  over  the  entire 
length  of  the  sample , In  Run  E-6,  this  may  have  been  m part  due  to 
the  fact  that  cutting  oil  was  used  during  the  final  sizing  cf  the 
PTFE  sample;  no  explanation  is  offered  in  the  case  of  Run  E-9. 

Test  Runs  E-7  and  E-10  involved  the  use  cf  solid  PTFE  rod 
fabricated  from  Teflon  9B  resin  by  the  Technical  Service  Laboratory 
at  duPont,  Cutting  oil  was  used  in  the  machining  of  the  PTFE  sample 
for  Run  E-7,  and  the  reaction  was  net  initiated  by  the  firing  cf  the 
igniter.  Test  Run  E-10  was  conducted  successfully  at  a PTFE-ilthium 
mass  ratio  of  0.219,  the  measured  reaction  propagation  velocity 
being  49  m/s. 

Reaction  pressure  measurements  were  made  during  Run  E-3 
with  Teflon  8A  rod,  Run  E-8  with  commercial  rod,  and  Run  E-9  with 
Teflon  8 red.  Since  the  propagation  of  the  reaction  in  Run  E-9 
was  incomplete,  no  meaningful  conclusions  may  be  drawn  from  this 
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pressure  measurement.  In  light  of  the  fact  that  strain  gauges  are 
very  temperature  sensitive,  it  is  felt  that  the  presence  of  a hot 
lithium  leak  near  one  of  the  strain  gauges  during  Run  E-8  rendered 
this  pressure  measurement  also  meaningless.  The  maximum  pressure 
measured  during  Run  E-3  was  19.3  bar,  the  indication  of  this  pressure 
occurring  well  after  the  reaction  front  had  passed  the  location  of 
the  thermocouple  nearest  the  cap.  It  was  therefore  concluded  that, 

t 

under  the  experimental  conditions  present  in  this  run,  as  shown  in 
Table  3,  the  pressure  wave  does  not  precede  the  reaction  front  through 
the  sample. 

From  the  results  of  the  test  runs  conducted  with  6.35  mm 
diameter  solid  PTFE  rod  samples , it  is  seen  that  in  no  case  did  the 
propagation  of  the  reaction  front  occur  normally  in  runs  involving 
PTFE  samples  on  which  cutting  oil  had  been  used  during  machining-. 

It  is  therefore  concluded  that  the  presence  of  small  quantities  of 
residual  oil  on  the  surface  of  the  PTFE  sample  in  some  manner  effects 
the  reaction  propagation  process.  Neglecting  the  results  of  Run 
E-l,  in  which  the  large  clearance  between  the  lioium  and  PTFE 
samples  was  present,  the  successful  test  runs  with  commercial  PTFE 
rod  were  carried  out  at  an  average  PTFE-lithium  mass  ratio  of  0.222, 
and  yielded  an  average  measured  reaction  propagation  velocity  of  49 
m/s.  The  average  PTFE-lithium  mass  ratio  for  the  successful  test 
runs  with  PTFE  rod  samples  fabricated  from  Teflon  8A  resin  was  0.226, 
and  the  average  reaction  propagation  velocity  was  46  m/s,  while  for 
the  successful  run  conducted  with  the  Teflon  9B  rod  sample,  the  PTFE- 
lithium  mass  ratio  was  0.219  and  the  measured  average  reaction 
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propagation  velocity  was  49  m/s.  Neither  of  the  test  runs  conducted 
with  Teflon  8 rod  samples  were  successful. 

It  is  believed  that,  within  the  limits  of  experimental 
error,  the  results  of  the  runs  with  the  commercial  rod  samples  may 
not  be  distinguished  from  the  results  of  the  runs  with  the  Teflon  8A 
rod  samples,  or  those  of  the  run  with  the  Teflon  9B  rod  sample.  For 
this  reason,  the  experimental  results  for  the  test  runs  involving 
6.35  mm  diameter  solid  PTFE  rod  samples  may  be  summarized  by  the 
average  values  from  the  successful  runs,  which  are  a PTFE-lithium 
mass  ratio  of  0.223,  and  a reaction  propagation  velocity  if  47.8 
m/s. 

4.4  Test  Runs  with  PTFE  Samples  of  Other  Diameters 

In  addition  to  the  test  runs  involving  9.52  mr  diameter, 
and  6 . 35  mm  diameter  PTFE  samples,  test  runs  were  conducted  with  3.18 
mm  (1/8  inch)  diameter,  4.76  mm  (3/16  inch)  diameter,  and  7.94 
mm  (5/16  inch)  diameter  PTFE  samples.  These  runs,  in  which  commercial 
PTFE  rod  was  used,  are  summarized  in  Table  4.  The  measuring  junctions 
of  the  thermocouples  used  in  these  test  runs  were  located  at  the 
lithium-PTFE  interface,  the  preparation  of  the  samples  was  carried  out 
in  the  inert-atmosphere  glove  box,  and  the  PTFE  samples  were  machined 
to  a diameter  nominally  0-25  mm  /0-010  inch)  smaller  than  the  core 
rod  diameter  to  a'’oid  an  interference  fit.  The  lithium  surrounding 
the  LiF  at  the  end  of  the  sample  was  neglected  in  the  computation  of 
the  PTFE-litbium  mass  ratios,  and  the  pressure  le'-ei  within  the  sample 
tube  at  the  time  of  the  firing  of  the  igniter  in  these  runs  was 
approximately  10  ^ bar  The  reaction  propagation  velocities  were 


Table 


>s 

<T3 

tc  -h  tn 
u ^ 

a.  c a 

O *-H  ^ 

W QJ 

fe  > 


\u>  to 

iUj  CO  Of) 

H ^ s- 

0-  r 


1^  CO  ool 
«■*  K >-  I 


OJ  F 
MSB 
Cl  flj  ''■» 

H 'H 

o.  a 


1 

1 

m 

CT» 

1 

l 

(N» 

" 

r* 

PN 

n* 

r-* 

f-4 

O 

O 

• 

O 

O 

O 

O 

H 

T—^ 

O 

O 

O 

O 

wO 

m 

vO 

vO 

^0 

<t 

CM 

<N 

VO 

vC 

v£> 

0 

r- < 

oo 

<T> 

00 

CO 

in 

in. 

lT 

1 r 

00 

00 

j£> 

r*v 

m 

>3 

vT 

f-H 

1 

m 

1 

CVJ 

1 

a: 

1 

•x 

1 

1 

’ — 1 

87 


determined  from  the  analog  output  records  from  the  CEC  Oscillograph, 
data  obtained  from  the  storage  oscilloscope  being  used  for  confirmation 
when  possible. 

As  shown  in  Table  4,  Runs  H-l  and  H-2  were  conducted  with 
3.18  mm  diameter  PTFE  rod  samples  under  similar  experimental  conditions, 
the  PTFE-lithium  mass  ratio  for  both  runs  being  0.047.  In  tnese  runs, 
the  reaction  propagation  velocity  measured  over  that  portion  of  the 
sample  between  the  thermocouple  nearest  the  igniter  and  the 
thermocouple  at  the  midpoint  of  the  sample  was  considerably  higher 
than  that  measured  over  the  portion  between  the  thermocouple  at  the 
midpoint  and  the  thermocouple  nearest  the  cap.  The  reaction  zone  was 
thus  observed  to  be  declerating  as  it  propagated  through  the  sample, 
it  being  therefore  concluded  that  normal  reaction  propagation  is 
not  possible  at  PTFE-lithium  mass  ratios  of  this  magnitude.  The 
pressure  within  the  sample  tube  during  these  runs  was  also  measured. 

A maximum  pressure  of  3.8  bar  in  Run  H-l  occurred  well  after  the 
reaction  front  had  passed  the  location  of  the  thermocouple  nearest  the 
cap,  indicating  that  the  pressure  wave  does  not  precede  the  reaction 
front  through  the  sample  in  runs  under  these  experimental  conditions. 
Although  the  reaction  front  did  propagate  over  the  entire  length  of 
the  sample  in  Run  H-2,  no  pressure  rise  above  the  initial  condition 
of  approximately  10  ^ bar  was  observed. 

Test  runs  J-l  and  J-2  were  carried  out  with  4.76  mm  diameter 
PTFE  rod  samples  under  similar  experimental  conditions,  at  PTFE-lithium 
mass  ratios  of  0.117  and  0.113,  respectively  As  shown  in  Table  u,  the 
reaction  propagation  measured  in  Run  J-l  was  25  m/S,  and  that  in  Run 
J-2,  19  m/s.  One  of  the  leadwires  from  the  thermocouple  nearest  the 
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igniter  was  broken  during  the  preparation  for  Run  J-l,  and  a similar 
occurrence  rendered  the  thermocouple  nearest  the  cap  inoperable  in 
Run  J-2.  In  these  two  test  runs,  the  reaction  propagation  velocity 
was  therefore  determined  usin®  data  olr--Jned  from  the  thermocouple 
located  at  the  midpoint  of  the  sample  tube,  and  the  one  other 
functional  thermocouple.  The  pressure  within  the  sample  tube  was 
measured  during  Run  J-2,  the  maximum  pressure  of  14.8  bar  occurring 
after  che  reaction  front  had  passed  the  location  of  the  thermocouple 
nearest  the  cap.  The  pressure  wave  thus  did  not  propagate  through  the 
sample  at  a more  rapid  velocity  than  did  the  reaction  front  in  this 
run.  Although  the  reaction  front  propagated  over  the  entire  length 
of  the  sample  in  these  two  test  runs,  the  character  of  the  propagation 
appeared  weak,  as  indicated  on  the  analog  output  record  from  the  CEC 
Oscillograph  by  a rather  slow  rate  of  temperature  rise  caused  by  the 
cnival  of  the  reaction  front  at  the  locations  of  the  functional 
thermocouples.  In  light  of  this  observation,  it  is  concluded  that 
the  propagation  of  the  reaction  is  marginal  under  the  conditions  of 
Runs  J-l  and  J-2,  as  shown  in  Table  4,  and  that  the  results  of  these 
two  runs  are  of  questionable  reliability. 

The  PTFE  rod  sample  utilized  in  Run  K-l  was  of  7.94 
mm  diameter,  as  shown  in  Table  4.  At  a PTFE-lithium  mass  ratio  of 
0.360,  the  measured  reaction  propagation  velocity  was  76  m/s,  this 
result  having  been  determined  using  the  data  from  the  two  thermocouples 
located  nearest  the  ends  of  the  sample.  Due  to  the  rather  high 
PTFE-lithium  mass  ratio,  local  melting  of  the  sample  tube,  accompanied 
by  sizeable  lithiun  leaks,  occurred  during  this  run  The  thermocouple 
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extension  vires  were,  however,  not  damaged,  and  the  reaction 
propagation  velocity  measurement  was  judged  reliable 

4.5  Experimental-Theoretical  Comparisons 

In  order  to  compare  the  predictions  of  the  model  with  the 
experimental  results,  it  is  first  necessary  to  numerically  evaluate 
the  characteristic  time  and  the  characteristic  velocity,  as  used  in 
Equation  (2.27).  The  characteristic  time,  t is  defined,  as 
described  in  Chapter  II,  as  the  time  required  to  thermally  decompose 
90X  of  an  initial  amount  of  PTFE  at  the  stoichiometric  adiabatic 
flame  temperature  at  one  atmosphere  total  pressure.  The  PTFE 
thermal  decomposition  rate  law  is  given  by  the  following  relationship: 


dZ  * 
T 

dt 


•L  /RT 


- Ao  V6 


(2.6) 


For  the  situation  in  which  the  temperature  is  constant  and  equal  to 
the  stoichiometric  adiabatic  flame  temperature  at  standard  conditions, 
the  product  of  the  frequency  factor  and  the  exponential  term  m 
the  above  equation  is  also  a constant,  and  is  defined  as: 

-E  /RT 
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(4.1) 


The  substitution  cf  Equation  (4.1)  into  Equation  (2.6)  yields: 


dZ£ 

dt 


K Z * 
T 


(4  2) 


Rearranging  and  integrating  the  above  equation,  and  applying  the 
initial  condition: 


t - 0 


Z * - 
T 


(4.3) 
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results  In  the  following  expression  for  the  undecomposed  PTFE  mass 
fraction  as  a function  of  time: 


£n  ZT*  - - Kt 


(4.4) 


Under  the  condition  that  the  thermal  decomposition  of  the  PTFE  is 
90%  complete,  the  value  of  Is  0.10,  and  the  corresponding  time 
is  the  characteristic  time,  as  given  by: 


_ - An(Q.lO) 
Cch  “ K 


(4.5) 


From  Reference  [37]: 

19  -1 

A - 3.0  x 10  s 
o 


and 

E - 3.475  x 103  J/mol 
a 

The  stoichiometric  adiabatic  flame  temperature  under  standard 
conditions  is  1966  t'  as  previously  described.  These  quantities  may 
be  combined  to  yield: 


E 

21.245  , (4.6) 

af 


and,  from  Equation  (4.1): 


K - 1.7805  x 10"10  s*1  . (4.7) 


The  value  of  t , is  then  determined  from  Equation  (4.5): 
ch 


t , - 1.293  x 10 
ch 
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It  is  now  desired  to  evaluate  the  characteristic  velocity  as  given 
by  the  following  expression: 


, -l/i 

ch  ” t 


T_ 

ch 


(2.27) 


Using  the  value  of  the  thermal  diffuslvity  of  PTFE  as  given  in 
Reference  [A3],  and  the  characteristic  time  as  evaluated  above,  this 
expression  is  evaluated  to  yield: 


V • 29.4  m/s 
ch 


(4.9) 


A summary  of  the  results  or  the  successful  llthium-PTFE 
test  runs  is  sho’/n  in  Table  5.  The  values  of  the  PTFE-lithium 
mass  ratios  and  the  measured  reaction  propagation  velocities  given 
in  Table  5 are  averages  for  all  successful  test  runs  with  PTFE 
samples  of  the  diameter  and  type  listed  with  the  exception  of  the  data 
for  the  7.94  mm  diameter  PTFE  sample,  with  which  only  one  successful 
test  run  was  conducted.  It  may  be  noted  from  this  table  that, 
neglecting  the  data  from  the  one  test  run  with  the  7.94  mm  diameter 
PTFE  rod  sample,  the  measured  reaction  propagation  velocities  are 
seen  to  decrease  with  increasing  PTFE-lithium  mass  ratic,  a trend 
which  was  predicted  by  the  model-  The  predicted  dimensionless 
reaction  propagation  velocities  shown  in  Table  5 were  evaluated  from 
the  linear  relationship  for  V^*  as  a function  of  PTFE-lithium  mass 
ratio,  Equation  (2.28),  anfl  the  dimensional  reaction  propagation 
velocities  were  determined  through  the  following  rearrangement  of 
Equation  (2.26): 
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The  error  values  given  in  Table  5 are  Che  absolute  value  of  che 
difference  between  che  measured  and  predicted  velocities  divided  by 
the  measured  velocity  for  each  of  the  PTFE  sample  diameters  and  types 
listed. 

In  the  numerical  solution  of  the  model,  the  properties  of 
the  PTFE  were  taken  as  those  of  the  solid,  as  previously  described. 

It  was,  therefore,  felt  that  the  experimental  results  from  test  run* 
in  which  solid  PTFE  rod  samples  were  used  should  provide  a more 
meaningful  test  of  the  model  than  che  results  of  the  runs  in  which 
PTFE  resin  samples  were  utilized.  From  Table  5 it  is  seen  that  for 
the  test  runs  involving  6.35  mm  diameter  PTFE  rod  samples,  the 
predicted  reaction  propagation  velocity  is  1.71  lower  than  the  average 
experimentally  measured  value.  In  contrast  to  this  close  agreement 
is  the  result  of  the  test  run  with  the  7.94  mm  diameter  PTFE  rod 
sample,  which  is  a measured  reaction  propagation  velocity  approximately 
441  higher  than  that  predicted  the  model.  No  explanation  for  this 
large  discrepancy  is  offered,  although  the  fact  that  only  one  test 
run  with  this  diameter  PTFE  rod  sample  was  conducted  casts  seme  doubt 
on  the  reliability  of  the  experimental  results. 

No  attempt  was  made  to  adjust  the  model  to  account  for 
situations  in  which  PTFE  resin  samples  were  used  since  the  estimation 
of  the  thermal  diffusivlty  of  a packed  granular  bed  yields  values 
which  are  approximate  at  best,  and  it  was  previously  found  that  the 
results  of  the  model  are  a strong  function  of  the  value  used  for  the 
thermal  diffusivlty  of  the  PTFE.  The  predicted  reaction  propagation 
velocities  for  the  test  runs  with  PTFE  resin  samples  were  arrived 
at  through  the  use  of  the  properties  of  solid  PTFE  in  the  model,  and 


were  included  in  Table  S for  comparison  purposes  only.  The  accuracy 
with  which  the  reaction  propagation  velocity  was  predicted  for  the 
test  runs  involving  PTFE  resin  samples  is,  therefore,  rather  surprising, 
in  light  of  the  uncertainty  in  the  value  of  o^,  as  well  as  the  fact 
chat,  for  the  6.35  mm  diameter  PTFE  resin  test  runs,  the  predicted 
value  of  V,,*  was  obtained  by  extrapolating  the  linear  relationship 
given  by  Equation  (2.28)  beyond  the  range  of  analytical  investigation. 

Under  experimental  conditions,  the  aspect  ratio  of  the 
lithium-PTFE  reaction  samples,  defined  as  the  length  of  the  sample 
divided  by  its  radius,  was  approximately  15.  From  the  results  of  the 
numerical  solution  of  the  model  previously  described,  however,  it  was 
indicated  that  the  predicted  dimensionless  steady  reaction  propagation 
velocity  was  not  a function  of  aspect  ratio  for  values  of  L*/H*  of  10 
and  greater.  It  is,  therefore,  believed  that  this  lack  of  geometric 
similarity  Is  not  a significant  source  of  error.  From  Table  5 it  may 
be  seen  that  the  reaction  propagation  velocity  was  predicted  mc-c 
accurately  for  the  test  runs  with  6.35  mm  diameter  PTFE  samples  than 
for  the  runs  involving  PTFE  samples  of  other  diameters.  A possible 
explanation  is  the  fact  that  the  model  was  cast  in  cartesian  coordinates, 
the  surface  to  volume  ratio  cf  the  PTFE  semple  being  therefore  only  a 
function  of  the  thickness  cf  the  PTFE,  which  was  held  constant. 

The  surface  to  volume  ratio  of  the  cylindrical  PTFE  samples  used  in 
the  test  runs,  however,  varies  inversely  with  diameter.  Additionally, 
it  is  estimated  that  an  error  in  the  axial  positioning  of  the 
thermocouples  of  approximately  + 3 mm  was  possible  during  the 
preparation  for  the  test  runs.  For  the  nominal  102  mm  instrumented 
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length  of  the  sample  tube,  this  positioning  error  would  bo  reflected 
as  a + 6%  error  in  the  reaction  propagation  velocity  measurements. 

The  estimated  total  error  In  these  measurements,  due  to  thermocouple 
positioning  errors  and  the  limitations  of  the  Instrumentation,  was 
+ 15  percent. 

The  purpose  of  the  pressure  measurements  made  during  several 
of  the  test  runs  was  not  to  precisely  determine  the  pressure  level 
at  which  the  reaction  took  place,  as  this  was  beyond  the  capability 
of  the  experimental  apparatus-  Rather,  these  measurements  were 
undertaken  in  order  to  assess  whether  the  reactants  ahead  of  the 
reaction  front  are  preheated  by  thermal  conduction  only,  as  postulated 
in  che  mechanism,  or  by  the  presence  cf  a hot  gas  pressure  wave  which 
precedes  the  reaction  front  through  the  sample.  The  bulk  of  these 
pressure  measurements  indicated  that  the  rejction  front  reached  the 
end  of  the  sample  opposite  the  igniter  prior  to,  or  coincident  with, 
the  arrival  of  the  pressure  wave,  lending  credence  to  the  assumed 
mechanism  of  heat  transfer  by  conduction  only. 

From  chese  comparisons,  it  is  therefore  concluded  that  the 
two-dimensional  model  based  on  the  reaction  mechanism  previously 
proposed  enables  predictions  of  the  reaction  propagation  velocity 
in  heterogeneous  lithlum-PTFE  samples  which  are  in  reasonable 
agreement  with  experimental  results.  To  argue  that  the  proposed 
reaction  mechanism  is,  in  fact,  an  accurate  rep  esentation  of  the 
situation  under  which  the  chemical  reaction  actually  takes  place 
would  be  Ludicrous.  It  is,  however,  observed  that  acceptably  accurate 
reaction  propagation  velocities  may  be  predicted  from  a model  based  on 


this  mechanism 


< 


CHAPTER  V 
SUMMARY 

The  overall  objective  of  this  study  was  the  investigation 
of  the  chemical  reaction  propagation  process  occurring  in  heterogeneous 
lithium-PTFE  reaction  samples.  This  process  was  examined  both 
analytically  and  experimentally  in  order  that  the  reaction  propagation 
velocity  might  be  predicteu  and  measured. 

A reaction  mechanism  was  proposed  in  which  the  reaction  rate 
controlling  step  was  assumed  to  be  the  thermal  decomposition  of  the 
PTFE  to  C^F^  gas  according  to  published  vacuum  pyrolysis  rate 
information.  It  was  postulated  that  the  C.,F^  produced  by  the  thermal 
decomposition  of  the  PTFE  reacts  immediately  upcn  its  formation  with 
the  lithium  to  form  condensed-phase  reaction  products  Conduction  of 
heat  from  the  reaction  zone  through  the  unreacted  lithium  and  PTFE 
was  assumed  to  be  the  sole  mode  ct  ene'gy  transfer  by  which  these 
reactants  are  preheated. 

Based  on  this  proposed  reaction  mechanism,  a two-dimensional 
model  was  developed.  This  mode.* , in  dimensionless  form,  was  solved 
numerically  tc  predict  the  steady  dimensionless  reaction  propagation 
velocity  as  a function  of  the  PTFE-lithium  mass  ratio.  Over  the 
range  analytically  inves tigaced , the  predicted  value  of  the  steady 
dimensionless  reaction  propagation  velocity,  V^* , was  observed  to 
decrease  with  increasing  PTFE-lithium  mass  ratio,  M* , according  to  the 
linear  relationship: 

V’s*  » 1.8575  - 1 1536  M* 

% 


\ 
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The  experimental  phase  of  this  study  involved  the  use  of 
cylindrical  reaction  samples  through  which  the  reaction  front  was 
allowed  to  propagate  axially.  The  cylindrical  PTFE  sample  was 
located  along  the  longitudinal  centerline  of  the  reaction  sample,  and 
«as  surrounded  by  the  lithium.  A test  apparatus  was  developed  which 
allowed  the  motion  of  the  reaction  front  to  be  monitored  by  means  of 
thermocouples  positioned  at  various  axial  locations  along  the  lithium- 
PTFE  interface,  the  reaction  having  been  initiated  at  one  end  of  the 
reaction  sample  by  the  action  of  an  electric  igniter  Reliable 
reaction  propagation  data  were  obtained  from  test  runs  in  which  PTFE 
samples  in  the  form  of  solid  rod  as  well  as  resin  were  utilized  at 
two  different  PTFE  sample  diameters.  The  propagation  of  the  reaction 
front  was  shown  to  be  axisymmetric  and  crude  pressure  measurements 
were  made  from  which  it  was  indicated  that  the  reactants  are  not 
preheated  by  gas-phase  convection.  The  estimated  accuracy  oi  the 
reaction  propagation  velocity  measurements  was  + 15  percent. 

From  comparison  with  experimental  results,  it  was 
demonstrated  that  the  model  predicts  reaction  prcpagation  velocities 
with  acceptable  accuracy  for  the  reaction  of  lithium  with  PTFE  in  the 
form  of  both  soiid  rod  and  ~usln  Better  correlation  was,  however, 
obtained  at  the  smaller  FTFE  sample  diameter. 

The  results  of  this  study  indicate  that  several  areas  are 
of  interest  for  future  investigation  In  order  to  obtain  more 
reliable  reaction  propagation  velocity  data  with  this  1 uel-oxidizer 
combination,  a mere  ilexible  high  speed  data  acquisition  system  with 
better  time  resolution  is  necessary.  Also,  an  igniter  capable  of 
consistently  initiating  the  chemical  reaction  under  ail  experimental 
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conditions  would  be  of  benefit.  The  range  of  PTFE-lithium  mass  ratios 
examined  with  the  model  should  be  extended,  especially  in  the 
direction  of  lower  values  of  this  parameter.  The  effect  on  the 
reaction  propagation  velocity  of  the  presence  of  an  additional  gaseous 
oxidizer  should  also  be  examined.  Finally,  experimental  investigations 
with  other  fuels  and  oxidi  ers  would  be  of  great  utility  in  determining 
the  general  validity  of  the  proposed  reaction  mechanism  and  the 
resulting  model. 
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APPENDIX  A 


DERIVATION  OF  THE  FINITE  DIFFERENCE  EXPRESSIONS 


During  the  transformation  of  the  dimensionless  heat  conduction 
equations.  Equations  (2  16)  and  (2.18),  into  finite  difference  form 
for  numerical  solution,  three  different  types  of  grid  patterns  were 
encountered.  These  three  patterns  are  shown  in  Figure  IS,  and  consist 
of  regular  grid  patterns  at  interior  grid  points.  Figure  15(a), 
grid  patterns  at  spatial  boundaries.  Figure  15(b),  and  grid  patterns 
with  unequal  grid  spacing  such  as  present  at  the  lithium-PTFE 
interface.  Figure  15(c),  The  derivations  of  the  finite  difference 
expressions  approximating  the  partial  differential  equations  for 
each  of  these  cases  are  illustrated  in  the  to) lowing  examples.  For 
ease  of  notation,  the  * superscripts  present  in  the  original  equations 
are  omitted;  superscripts  in  the  following  expressions  refer  to 
dimensionless  time  level  and  subscripts  to  spatial  location 

For  regular  g*-id  patterns  at  interior  grid  prints  in  the 
stcichiocetric  lithium,  as  shown  in  Figure  15(a),  finite  difterence 
expressions  approximating  che  solution  or  Equation  (2,18), 
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(2 . 18) 


may  be  arrived  at  by  first  writing  a Taylor's  Series  expansxon  about 
the  grid  point,  i,  j,  as  follows,  all  derivatives  being  evaluated  at 
grid  point  i , j : 
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It  Is  desired  to  add  the  above  four  equations  In  such  a manner  as  to 
arrive  at  a finite  difference  expression  approximating  the  spatial 
derivatives  on  the  right-hand  side  of  Equation  (2.i8).  Accordingly, 
Equations  (A.l)  through  (A. 4)  are  multiplied  by  the  constants 
Cf,  and  , respectively.  The  following  constraints  are  then 

necessary  in  order  for  the  addition  to  yield  the  desired  results: 


^coefficients  of 


a t 
ax2 
3T 


£ coefficients  of  • 0 , 

92T 

^coefficients  of  ~~y  “ 1 , 

3Y 

Y coefficients  of  ~ • 0 


(A  5) 


These  constraints  applied  to  the  addition  c:  Equations  (A  1)  through 
(A  4)  then  yield: 


fr  <ci  + V ‘ 1 • 


AX(CX  - C2’>  - 0 , 
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— <c3  4 c4)  - 1 , 


AY,  fC,  - C,  > • 0 
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(A. 6) 
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Solution  of  Equation  (A. 6)  yields: 


c!  - c2 ^2 

1 * AX* 


c3  “ c3  " ~T 

J J A Y * 


(A.;) 


Multiplying  Equati:ne  (A.l)  through  (A. 4)  by  the  constants  through 
as  given  by  Equations  (A. 7),  respectively,  adding  the  resulting 
expressions,  and  neglecting  terms  involving  the  fourth-and-higher 
order  derivatives  gives: 

4*24- jJt*  , * r , -n"  . _Ut»  ,.1  + t » , 2i^  I. 


(A  . 8) 


The  dimensionless  time  derivative  on  the  left-hand  side  of  Equation 
(2.18)  is  then  approximated  by  the  forward  difference  expression 


if;  - tn  , 
..  JL?J. 


(A. 9) 


Multiplying  Equation  (A- 8)  by  the  dimensionless  coem.lent  A*,  and 
equating  the  resulting  expression  with  Equation  (A  9)  yields  the 
following  f inice  difference  approximaticn  for  Equation  (2  18): 


-^ilTN  +T* 

AX2  [Vl.j  Ti-l,j 


2XN  +TN  -2TN  . 

i.jj  AY  2^i.j«t-l  ‘i.J-l  ,j|’ 

(A. 10) 


This  equation  may  then  be  rearranged  to  give  an  explicit  expression 
for  the  dimensionless  temperature  at  grid  point  i,j  at  dimensionless 


F 
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time  level  N+l,  in  terms  of  the  dimensionless  temperatures  at  grid 
point  i,j  and  the  four  adjacent  grid  points  at  dimensionless  time 
level  N a^  follows: 


tN*‘  . 

AX 


i.J 


....  . tN  ,tN 

fTi+i.j  Ti-i,j ' 2Tij 


♦ J**!.!  • <*•“> 


The  finite  difference  approximation  for  Equation  (2.16),  which 
includes  a heat  generation  term,  at  interior  grid  point  i,j  In  the 
PTFE  is  derived  in  a like  manner  to  give: 

A.  , At 


- ti*l _1  tn  ♦ tn  - 2TN 
i.J  2 (Vl.j  ‘i-l.j 


dZ 
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AX 
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(A. 13) 


For  grid  points  on  spatial  boundaries  as  shown  in  Figure 
15(b)  which  are  assumed  to  be  adiabatic,  finite  difference  expressions 
approximating  the  solution  of  Equation  (2  18)  are  derived  by  again 
writing  a Taylor's  Series  expansion  around  grid  point  1 , j in  the 
stoichiometric  lithium  as  follows,  ail  derivatives  being  evaluated 
at  grid  point  i , j : 


TN  - tn 
i+l.J  i.j 


2 2 3 3 

3T  . AX^  31  , AX  31 


(A. 14) 
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(A. 16) 


Noting  chat  at  an  adiabatic  boundary  as  shown  In  Flgura  15(b), 

3T 

• 0,  and  naglacting  terms  involving  thitd-and-hlgher  order 

derivatives  allows  Equation  (A. 14)  to  be  rearranged  tc  yield  an 
32T 

expression  for  — r.  It  being  desired  tc  add  Equations  (A. 15)  and 
3X  32T 

(A. 16)  to  give  an  expression  for  — =■,  these  equations  are  first 

multiplied  by  constants  and  C^,  respectively  The  vali-es  of 

these  constants  are  determined,  through  constraints  on  the  sums  of 

2 

the  coefficients  of  — and  r^,  to  be: 

3Y2  3Y 


C.  - C - — ~ 
5 6 ayl2 
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Multiplying  Equations  (A. 15)  and  (A. 16)  by  and  as  given  above, 

adding  the  resulting  expressions  while  neglecting  f cut th-and-higher 

d2T 

oraer  derivatives,  and  introducing  the  expression  for  — r yields: 

3X 

32T  . 92T  2 |_N  _N  \ ._1_/tN  . -N  ,.N  \ 

dX2  3Y2  " AX2  l 1+1  ' ‘'Jj  AY  2 I i*J+1  ".S*1  ' i*J/  ’ 


(A. 18) 


for  a grid  point  i,j  on  a spatial  boundary  as  shewn  in  Figure  15(b). 
Expressing  the  dimensionless  time  derivative  as  given  in  Equation  (A. 9), 
multiplying  Equation  (A. 18)  by  the  dimensionless  parameter  A*,  equating 
these  two  expressions  and  rearranging  gives  the  following  equation 
for  the  dimensionless  temperature  of  boundary  grid  point  i,j  in  the 
stoichiometric  litnium  at  the  N+l  time  level: 
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In  a similar  manner,  a finite  difference  approximation  for  Equation 
(2,16)  at  a spatial  boundary  grid  point  i,j  in  the  PTFE  layer  may  be 
derived  to  give: 
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For  grid  points  in  the  stoichiometric  lithium  which  are 
adjacent  to  the  lithium-PTFE  interface,  as  shown  in  Figure  15(c), 
finite  difference  expressions  approximating  the  solution  of  Equation 
(2.18)  may  also  be  arrived  at  by  writing  a Taylor's  Series  expansion 
about  grid  point  i,j  as  follows,  all  derivatives  being  evaluated  at 
grid  point  i ,j : 
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It  being  again  desired  to  add  the  above  four  equations  to  arrive  at  a 
finite  difference  approximation  for  the  spatial  derivatives  on  the 
right-hand  side  of  Equation  (2.18),  Equations  (A. 21)  through  (A. 24) 
are  multiplied  by  constants  through  C^q,  respectively.  The  values 
of  these  constants  are  then  determined  by  thr  constraints  given  by 
Equations  (A. 5)  to  be: 


C8-‘-J  * 


ayltayi(1,  + sv^r)  * 


Multiplying  Equations  (A. 21)  through  (A. 24)  by  constants  C7  through 
C^y,  respectively , adding  the  resulting  expressions,  and  neglecting 
terms  involving  thlrd-aud- higher  oiler  derivatives  yields: 
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where  the  values  of  and  C.^  are  given  by  Equations  (A. 25) 
Multiplying  Equation  (A. 26)  by  the  dimensionless  parameter  A*, 


Ill 


representing  the  dimensionless  time  derivative  In  Equation  (2.18)  as 
shown  In  Equation  (A. 9),  equating  and  rearranging  gives  the  following 
expression  for  the  dimensionless  temperature  at  grid  point  i,j  in  the 
stoichiometric  lithium  adjacent  to  the  lithlum-PTFE  Interface  at 
dimensionless  time  level  N+l: 

ax' 
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where  and  are  given  by  Equations  (A. 25).  Similarly,  a finite 
difference  approximation  for  Equation  (2.16)  at  grid  point  i,j  In  the 
PTFE  adjacent  to  the  llthium-PTFE  interface  may  be  derived  to  yield: 
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The  equations  used  in  the  numerical  solution  of  the 
dimensionless  form  of  the  two-dimensional  model,  Equations  (2.16) 
through  (2.20)  were  either  identical  to  those  presented  in  this 
appendix,  or  combinations  of  two  of  the  examples  illustrated. 
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APPENDIX  C 


PHYSICAL  AND  CHEMICAL  PROPERTIES 

The  values  of  the  physical  and  chemical  properties  used  in 
the  numerical  solution  of  the  model  are  summarized  in  Table  6. 


Table  6 


Values  of  Physical  and  Chemical  Properties  Used  in  the  Model 


Property Value Source 


A 

o 

3.0  x 1019  s"1 

Reference  [37] 

CP 

4.187  x 103  J/kg-K 

Reference  [42] 

*L 

% 

2.470  x 103  J/kg-K 

Reference  [1] 

% 

1.047  x 103  J/kg-K 

Reference  [43] 

E 

a 

3.475  x 105  J/mol 

Reference  [37] 

ahr 

1.237  x 10^  J/kg  reactants  Computed  using  References 

and  [37) 

[1] 

55.45  tym-K 

Reference  [42] 

*P 

1.731  W/m-K 

Reference  [44] 

kT 

0.2527  W'rn-K 

Reference  [43] 

Taf 

1966  K 

Computed  using  References 
[2]  and  137] 

tl), 

W 

1.416 

Conputed  using  Reference 

[1] 

PL 

481  kg/m3 

Reference  [42] 

PP 

1110  kg/ra3 

Computed  using  References 
[42]  and  [43] 

[1], 

P j 

2160  kg/m3 

Reference  [43] 
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